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A theory of intense-field dynamic alignment and high harmonic generation from 
coherently rotating molecules and interpretation of intense-field ultrafast pump-probe 

experiments 

A. Abdurrouf 1 -^ and and F.H.M. Faisal^ 

1 Fakultat fiir Physik, Universitdt Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany and 
2 Department of Physics, University of Brawijaya, Malang, Indonesia 65145 

A theory of ultra-fast pump-probe experiments proposed by us earlier [F.H.M. Faisal et at, Phys. 
Rev. Lett. 98, 143001 (2007) and F.H.M. Faisal and A. Abdurrouf, Phys. Rev. Lett. 100, 123005 
(2008)] is developed here fully and applied to investigate the phenomena of dynamic alignment and 
high harmonic generation (HHG) from coherently rotating linear molecules. The theory provides 
essentially analytical results for the signals that allow us to investigate the simultaneous dependence 
of the HHG signals on the two externally available control parameters, namely, the relative angle 
between the polarizations, and the delay-time between the two pulses. It is applied to investigate 
the characteristics of high harmonic emission from nitrogen and oxygen molecules that have been 
observed experimentally in a number of laboratories. The results obtained both in the time-domain 
and in the frequency-domain are compared diwith the observed characteristics as well as directly 
with the data and are found to agree remarkably well. In addition we have predicted the existence 
of a "magic" polarization angle at which all modulations of the harmonic emission from nitrogen 
molecule changes to a steady emission at the harmonic frequency. Among other things we have 
also shown a correlation between the existence of the "magic" or critical polarization angles and 
the symmetry of the active molecular orbitals, that is deemed to be useful in connection with the 
"inverse problem" of molecular imaging from the HHG data. 
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PACS numbers: 32.80.Rm, 32.80.Fb, 34.50.Rk, 42.50.Hz 



I. INTRODUCTION 



In recent years there has been much interest and 
progress in understanding the interaction of atoms and 
molecules with intense laser fields (e.g. reviews [l|, Hj]). 
Among the phenomena observed, the high-order har- 
monic generation (HHG) is of particular interest, no 
less because of its potential applications as a source of 
coherent ultraviolet light and/or for generation of ul- 
trashort attosecond laser pulses. In contrast to atoms, 
molecules have extra degrees of freedom such as vibra- 
tion, and rotation of the molecular frame, and have addi- 
tional symmetry properties, that give rise to richer phys- 
ical phenomena when they interact with intense laser 
pulses. Among them is the phenomenon of alignment of 
linear molecules by strong and longlaser pulses which 
has been investigated in the past H, IE @- Much 
interest has recently been generated by the observa- 
tion of recurrent dynamic alignments of linear molecules 
like N2 and O2 @, interacting with intense ultra- 
short laser pulses. They are monitored, for example, by 
non-destructive high harmonic generation signals from 
intense-field pump-probe experiments with delayed pairs 
of intense ultrashort pulses @, ES El El El . The dy- 
namic HHG signals have been used also to "reconstruct" 
the molecular orbitals [3, EE EH , to investigate proton 
motions [ItJ and molecular dynamics [l8l ]. 
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In this paper we derive fully a recently proposed |19l . 
I20I ] quantum theory of intense-field dynamic alignment 
and high harmonic generation from linear molecules and 
apply it to analyze the observed dynamical HHG signals 
for N2 and O2 molecules. Theoretical expressions for the 
signals are given analytically as a simultaneous function 
of the two external operational parameters - the delay 
time, td, and the relative polarization an gle, a, between 
the pump and the probe pulse [1, EE EH EE EH • 

Before proceeding further, we briefly discuss the main 
experimental characteristics of dynamic alignment and 
the HHG signals as observed for N2 and O2. We recall 
at the out set that the quantum measure of dynamical 
alignment of a rotating molecule is the quantum expec- 
tation value (with respect to the rotational wave-packet 
states induced by the pump pulse) of the "alignment 
operator" cos 2 6, that is averaged over the Boltzmann 
distribution of the initially occupied rotational states: 
A(td) = ((cos 2 9)) (td), where 9 is the angle between 
the molecular axis and the probe polarization direction; 
the double angular brackets stand for the expectation 
value with respect to the rotational wave-packets (inner 
brackets) and the statistical average with respect to the 
Boltzmann distribution (outer brackets) of the initially 
occupied rotational states. It was observed experimen- 
tally [E EE EH, El, EH that the dynamic (or delay-time 
dependent) HHG signal for N2 mimicked the "alignment 
measure" A (td). It exhibited the phenomenon of rota- 
tional revivals [E [E [2l[ including the "full-revival" with 
a period T r = jhc^ wnere B 1S the rotational constant 
[22| . as well as a ^-revival, and a 4-revival. They are 
consistent with the time dependence of A (td) defined 
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above, since the operator cos 2 9 can couple the rota- 
tional states with A J = ±2 (Raman allowed transitions) 
among the rotational states of the induced wavepackets, 
and thus can give rise to fractional revival periods associ- 
ated with the corresponding beat frequencies. In the case 
of O2, unexpectedly, an additional ^-revival appeared in 
the HHG signal [ill, E3, El- The latter is impossible 
for the alignment measure A (td) to account for, since 
it can not couple the rotational states with AJ = ±4, 
that could give rise to a beat period |T r . Thus, to fit 
their data of O2, Itatani et al. ll] proposed, empiri- 
cally, to consider the expectation value of the operator 
B (td) = ((sin 2 29)) (td). Subsequently, some of the early 
theoretical models of the HHG signal (e.g. [H, 0, HI) 
gave a similar result for O2 and thus appeared to jus- 
tify the empirical fit. Such a model also suggests that 
the maximum HHG signal for N2 can occur when the 
field polarization and the molecular axis were parallel, 
whereas the maximum signal of O2 would occur when 
they are "diagonal" (i.e. make an angle 9 = 45°). Un- 
like the time dependent signals themselves, their Fourier 
transform (F.T.), with sharply defined individual spec- 
tral lines and series, provide an alternative (and rather 
more precise) means of studying the dynamic alignment 
phenomenon. More recent experimental observations of 
the dynamic HHG signals for N2 and O2, and their F.T. 
have revealed surprising characteristics that can not be 
fully understood in terms of the earlier considerations. 
Thus: 

(a) Kanai et al. [l2j found that their experimental HHG 
signals for N2 and O2 could not be well fitted, respec- 
tively, by the expectation values of the operators cos 2 9 
and sin 2 29, alone. They considered empirically addi- 
tional operators involving higher powers of cos 2 9, or Leg- 
endre polynomials, to fit their data. 

(b) Miyazaki et al. [l3| measured the dynamical HHG 
signals of N2 and O2 and Fourier transformed their sig- 
nals and found not only spectral series containing strong 
Raman allowed but also weak Raman forbidden and 
anomalous lines, for both N2 and 02- 

(c) Itatani et al. [3] observed that the HHG signal from 
dynamically aligned N2 was enhanced when the pump 
polarization was taken parallel to the probe polarization, 
and were suppressed when the polarizations were taken 
to be perpendicular. 

(d) Kanai et al. [l2j and Miyazaki et al. [H, [26[ mea- 
sured the HHG signal for the diatomic N2, O2, and the 
triatomic CO2, for different relative angles a between 
the pump and probe polarizations, and observed that 
the HHG signal modulations are not only smaller in the 
perpendicular case, compared to the parallel case, but 
also are of oppositevhase in the two geometries. 

(e) Kanai et al. \12\ proposed a planar emission model 
of HHG which produced an opposite phase relation, as 
observed, but it did not yield the unequal modulation 
amplitudes, observed in the two geometries. 

(f) The present theory predicted (cf. and below) a 
"magic" polarization angle, a c « 55°, at which the har- 



monic emission from coherently rotating molecules with 
a g orbital symmetry (e.g. N2) becomes equal for all de- 
lay times td- Most recent observations by Yoshii et al. 
[271 [5^ appear to confirm the same. 

In this paper we present an ah initio development of 
the above mentioned theory [lj| [2(| that is shown to 
provide a unified theoretical account of all the phenom- 
ena noted above and other related characteristics of dy- 
namic alignments and the HHG signals from the aligning 
molecules, as well as their Fourier spectra, that have been 
observed experimentally. To this end, below we begin 
with a short schematic description of a typical intense- 
field pump-probe experiment on dynamic alignments and 
the molecular HHG signals as a function of (a) the time 
delay, td, and (b) the relative polarization angle, a, be- 
tween the pump and the probe pulse. In Sec. II, III 
and IV, we systematically derive the S-matrix theory of 
molecular alignment and dynamic HHG signal from an 
ensemble of freely rotating linear molecule, discuss the 
connection between the "one" - and the "many" -molecule 
signals, and the relation between the quantum amplitude 
for the emission of the HHG photons and the expecta- 
tion value of the dipole transition moment. In Sec. V, 
we apply the theory to N2 and O2 molecules and ob- 
tain analytic expressions for the "HHG operators" and 
the HHG signals for an arbitrary, a, and delay times, td- 
In Sec. VI. we use the theoretical expressions to explic- 
itly calculate the HHG signals for specific experimental 
parameters, for the linear molecules N2 and O2, both 
in the time domain and in the frequency domain, and 
discuss the results with reference to the corresponding 
experimental observations. In addition, we investigate in 
section VII a number of related problems of general in- 
terest including the influence of the probe pulse on the 
dynamic alignment, the effect of the initial temperature 
on the HHG signal, and the mean energy of the molecule 
after interaction with the pump pulse. We also discuss 
two other definitions used earlier for the HHG signal, as 
well as investigate the case of "adiabatic alignment" of a 
linear molecule, within the present theory, in the limit of 
long pulse durations. We end with a concluding summary 
in Sec. VIII. 

In Fig. [T]we show a schematic of a typical intense-field 
pump-probe experiment. A laser beam is first split into 
two parts, Li and L2, by a beam splitter (BS) with a 
desired ratio of the beam intensities. The probe-pulse L2 
is delayed by passing through a delay line system (D), 
by a finite amount td, with respect to the pump-pulse L\ 
and both are sent through a beam mixer (BM) to the 
target gas molecules from a gas jet. The high harmonic 
signal produced by the probe pulse is recorded by the 
detector system for each selected values of td- In addi- 
tion, a polarizer P can be inserted to rotate the angle of 
polarization of the probe pulse with respect to the polar- 
ization direction of the pump pulse at any desired angle 
a. The pulses are generally assumed to be effectively 
non-overlapping (td ^ 0) and that they are shorter than 
the period of the rotational degrees of freedom of interest. 
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Note that both td and a provide controllable parameters 



H_ 




BM d F 



is the same as the initial state) with respect to the tar- 
get molecule, and (for an ideal gas medium) along the 
forward direction of the incident field (cf. @, §4], and 
references cited in that section). 

Thus, we may write the laser- molecule interaction 
Hamiltonians appearing above as given by 



V N - Ll (t) = -M-Fx(t) - -F x (^(i)) : a : F x {4> 3 {t)) (2) 

where fj, is the permanent dipole moment (if non-zero) 
and a with Cartesian components a^; («,«') = (1,2,3) 
is the polarizability tensor of the molecule (always non- 
zero); and 



Figure 1: A scheme of a typical pump-probe experiment. See 
text for further explanation. 



V e -L 2 (<f>j{t - t d )) = -d e F (&(f - t d )) (3) 
where d e stands for the electronic dipole operator. 



II. A QUANTUM THEORY OF INTENSE-FIELD 

PUMP-PROBE EXPERIMENTS AND 
MOLECULAR HIGH HARMONIC GENERATION 
SIGNALS 

A. Total Hamiltonian and Equations of Motions of 
the Dynamical System 

Within the adiabatic Born-Oppenheimer approxima- 
tion of the target molecule, the total Hamiltonian of the 
system can be written [19j : (in a.u.: e = h = m — ac — 
1) as, 

H t ot{t) =H N + V N - Ll {t) +H e + V e - L2 {t - t d ) (1) 

where Hn is the nuclear Hamiltonian, VN-L t {t) is 
the interaction due to the pump pulse with the nu- 
clear motion at time t, H e is the electronic Hamilto- 
nian, and V e -L 2 (t ~~ td) i s the interaction of the probe 
pulse with the active electron at a delay td. Wc 
describe the two laser pulses (in the long-wavelength 
dipole-approximation) of the laser fields, F(<pj(t)) = 
i(t) cos (4>j(t)) and the corresponding the vector poten- 
tials by A((j)j(t)) = -^f(t) sin ((/>,■(*)), where f(t) is the 
slowly varying envelope of the electric field (compared to 
the period of the high harmonics, or the electron motion). 
The phase 4>j (t) of the field at the position of the active 
electron of the molecule is given by <j)j{t) = (uit — k^.Xj), 
where the CM. of the molecule is assumed to be located 
at a position "Xj " ; to and k w are the laser of frequency 
and the wavenumber, respectively. For the sake of sim- 
plicity of writing, we may suppress the notation the full 
Xj and t dependence of the phase factor <f>j (t) unless oth- 
erwise needed explicitly, e.g. while summing coherently 
the " many-molecule" emission amplitudes from different 
locations {Xj} to obtain the total amplitude associated 
with the signal macroscopically coherent signal, ft will be 
found that the coherent signal appears significantly for 
the "elastic" scattering (the final state of the molecule 



B. Total Wavefunction in Intense-field S-Matrix 
Theory 

We first consider a systematic solution of the time- 
dependent Schrodinger equation of the system 



iJU(i) = H tot (t)*(t) 



(4) 



using the general technique of intense-field many-body 
S- matrix theory (IMST) @. In this approach the total 
wavefunction of the system satisfying a given initial (fi- 
nal) condition can be written as a series expansion in such 
a way that the dominant virtual states, when present, 
can appear already in the leading terms of the series. To 
this end we introduce three partitions of the same to- 
tal Hamiltonian, referring to the initial, "i", the final, 
"/" , and the (deemed to be relevant) intermediate vir- 
tual state, "0" , interactions plus the corresponding "ref- 
erence" Hamiltonians: 



H tot {t) 



Hi + Viit) 
H f (t) + V f (t) 
H (t)+V (t). 



(5) 



It is also useful to define the reference Green's functions 
associated with the reference Hamiltonians, H s (t); s = 
i,f,0 : 



d_ 

dt 



H s (t))G s (t,t') = 5(t-t'). 



(6) 



In general, the Green functions can be obtained from 

' Mi 



the complete set of the fundamental solutions, ip) v) 

of the Schrodinger equations governed by the reference 
Hamiltonians H s (t); s = i, /, 0: 

G s (t,t') = -ie(t-t')J2\4 s) ( t )){^ s) ( t ')\ ( ? ) 

all j 
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The validity of the solutions Eq. ([6]) can be readily es- 
tablished by operating on the left hand side of Eq. 
with {ijft — Ho(t)), using Eq. ((6]) and the completeness 

of the fundamental solutions, J2j yPj S \t) = 1 

and the relation §^0{t — t') = S(t — t'), to obtain a delta- 
function integration on the right hand side, followed by 
the obvious simplification. Thus, we can express the to- 
tal wavefunction of the interacting system, evolving from 
an arbitrary initial state, |Xi(i)), as a series: 



with 



3=0 



*i 0) (*)) = \Xi(t)) 



(8) 



(9) 



*j X) (*)) - V, (tx) | Xi (*i)) (10) 



(a) i 



(b) i 



{k} + G 

X 



{M 



Figure 2: Quantum amplitude for coherent emission of a high 
harmonic photon (frequency fi) is the sum of a direct (a) and 
a time-reversed (b) diagram. 



same result for the single photon HHG emission ampli- 
tude can also be obtained using the ordinary quantum 
mechanics, simply by taking the quantum electrodynam- 
ically normalized interaction V*(i) for the spontaneous 
emission of a photon of frequency fl and wavevector 
(cf. Lecture 2]): 



7*(i)=JVne i *<We n -d e 



(13) 



*f (*)) = 



and 



dt 2 dhG u f (t,t 2 )V f {t 2 )G (t 2 ,h) 
U Jti 

xVi(tx)\xi(h)) 

(ii) 



dt n dhdhGj (t,t n ) 

xV f (t n ) Go (t 3 ,t 2 )V f (t 2 )G (t 2> h) 

xViMlxifa)). (12) 



In the above, Nn = 



T 3 " 



III. MANY-MOLECULE VS. 

SIGNALS 



ONE-MOLECULE 



Transition Amplitudes for High Harmonic 
Generation 



Emission of a harmonic photon of frequency fl = ncu 
and wavevector Kq, from its vacuum state |0n) (zero oc- 
cupation number in Fock-space), into a singly occupied 
number state, |ln), is fundamentally a quantum electro- 
dynamical process i.e. due to the interaction of the active 
electron with the vacuum-field albeit in the presence of 
the intense external laser field. Its theoretical formula- 
tion therefore clearly requires one to consider at least the 
combined state of the interacting "laser field (semiclassi- 
cal) + molecule + vacuum- field"- system in the extended 
space consisting of the direct product of the ordinary 
space of "laser field (semiclassical) + molecule" and the 
occupation number space of the vacuum and the emit- 
ted photon (cf. e.g. 0, §4.5]). Nevertheless, exactly the 



I? is the quantization vol- 
ume, Cfj is the polarization vector of the emitted photon, 
and d e is the usual electronic transition dipole opera- 
tor; the phase $>j(t) = (fit — k u ■ Xj). As usual in the 
present dipole approximation, we have neglected the re- 
tardation factor, e~ iKa ' r Rj 1; we may note explicitly 
that the exact position of the electron with respect to an 
arbitrary coordinate origin is given by Xj + r, where as 
before Xj is the CM. of the jth molecule and r is the 
position of the electron with respect to the CM. of the 
molecule. 

The HHG amplitude for the emission of a harmonic 
frequency f2 from the jth molecule is given by (cf. 0, 
§4]) by the sum of two "Feynman-like" diagrams, (a) 
and (b), shown in Fig. O The diagram (a) corresponds 
to the so-called "direct" amplitude (associated with the 
retarded Green's function), whereas the diagram (b) cor- 
responds with the "time-reversed" amplitude (associated 
with the advanced Green's function). The amplitude for 
the harmonic emission process can be written down an- 
alytically from diagram (a) by reading in the forward 
(upward) direction of time: First, the molecule is pre- 
pared by the "pump" laser pulse in the state i = |$,(t)). 
Next, the molecule interacts with the probe laser field 
(short horizontal line with a cross), then it propagates 
through Go = Gq (vertical line) . Next it interacts with 
the vacuum-field by yM (horizontal line ending in fl), 
and emits the harmonic photon of frequency fl, and fi- 
nally returns to the same state i = |<frj(t)) as before. 
A similar interpretation holds for the time reversed dia- 
gram (b) except that here the system evolves backward 
in time thorough Go = G^~\ The quantum amplitude 
Aj(fl), for the emission of the HHG photon of frequency 
f2, and wavenumbcr from the molecule located at 
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Xj, is given by the sum of the two diagrams: 

Aj(ST) = diag.(a) + diag.(b) (14) 
Mathematically, we have 

/OO />OC 
dt / dt' (4>i(t)\ V*(t) 
-oo J — OO 



xG [ +> (t,t')V e - L2 (t' ~t d )\ Xl (t')) (15) 



and, 



OO f oo 



diag.(b) = -i dt dt' (^i(t')\V e - L2 (f - t d ) 

J — cc J —oo 



x<V(i',t)F*(i)| X< (*)> 



(16) 



B. A Relation between HHG Amplitude and F.T. 
of Dipole Expectation Value 

i. Recasting the quantum HHG amplitude 

Before proceeding further, we consider the relation be- 
tween the quantum HHG amplitude, Eq. (| 14[) . and the 
expectation value of the electric dipole operator, d e , that 
is popularly used for calculations of HHG signals. To this 
end we first rewrite the quantum amplitude Eq. (|15[) by 
introducing the first order wavefunction (cf. Eq. 

(HHJ)) that arises from the initial state \xi(t)) due to the 
interaction with the probe pulse: 



= / dt'G^\t,t') 

' J —oo 

xV e - L2 (t ! -t d )\ Xi (t')) 
We may first rewrite Eq. (fl"5"|) as: 

diag.(a) = -i [ dt ( Xi (t) \V* (t)\ 
dte i{nt-K a . Xj ) 



(17) 



Nne n ■ d t 



(18) 



where we have used the explicit form of the interaction 
V*(t). Next we rewrite Eq. p^]) using a standard re- 
lation satisfied by the advanced and the retarded Green 
functions (e.g. [23|): 



Gi-\t',t)=\Gi + \t,t>) 



(19) 



We also note that the laser-molecule interaction is real 
(Hermitian), i.e. 



Ve-L 2 {t) = [Ve-L 2 {t)T 



(20) 



Thus, the integral over dt' in Eq. (| 16[) can be rewritten 

as, 



diag.(b) 



dt 



dt' {Xi(t')\V e - L2 (t'-t d ) 



xG^\f,t))V*(t)\ Xi (t)) 



-i I dt 

-oo 



dt' 



x 'G^\t,t')V e -L 2 (t'-t d )\ X i(t')) 
xV*(t)\ Xi (t)) 

dt(*M(t)\V*(t)\ X i(t) 



Xi(t) 



dte i(nt-K n . Xj 



(21) 



Hence, adding Eqs. (JTHJ) and (|2"Tj) , we get the quantum 
HHG amplitude in the suggestive form: 



Aj(Q) = -Waen- / dte^ nt ~ K ^ x ^ 

J — OO 

x {(xi(t) \d e \ + \d e \ Xi (t))} 

/OO 
-oo 

x{( Xi (i)|d e |*«(i))+c. C .} (22) 
where "cc." stands for the complex conjugate. 

2. Dipole expectation value 

The expectation value of the dipole operator, Dj $(t), 
of the transition dipole operator can be calculated within 
the lowest order KFR approximation [3(J HH, H2] of the 
wavefunction of the system as follows: 



Di,i(t) = (V(t) d e 

= (xi(t) de Xi(t) 



d, 



m)} 



= 1 4 1 +c.c.}+--- (23) 

Note that the zeroth order term in the first line above 
vanishes for centrosymmetric systems; we have also ne- 
glected the quadratic powers of the first order KFR- 
correction (and the higher order terms). 

Combining Eq. (J53J with Eq. ![22"]). the quantum HHG 
amplitude can be expressed in the form 



Am 



-i e 



-iKn-Xj 



/oo 
dte mt D iA (t) (24) 
-OO 
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Thus, the quantum HHG amplitude is clearly propor- 
tional to the Fourier transform (F.T.) of the expectation 
value of the transition dipole operator, Eq. ([23jl . We 
may recall that the proportionality constant Nq above 
is of quantum electrodynamical origin and can not be 
derived from the classical electrodynamics alone [33| . 

In practice, the F.T. of interest can be conveniently ob- 
tained by Fast Fourier transform (or FFT) numerically 
[3~4j . Alternatively, for "slowly varying" pulse envelopes 
(compared to the high harmonic frequency) one may ex- 
press the F.T. of Di t i(t) as a Fourier series (35[: 



Di,i{t) 



-in(wt— fe w -Xj 



D{nuj) 



(25) 



where D(nu) is the F.T. evaluated at the nth harmonic 
frequency f2 = nu. Thus, in terms of the F.T. compo- 
nents, the HHG amplitude Aj(Q) becomes, 



Aj(Q) = V-27ri£(0-nu/)e 



xN n en ■ D iti (nu) 



i(Ka—nk^)-Xj 



(26) 



where we have carried out the time integration over dt, 
in terms of the Dirac delta-function. 



C. Coherent Sum of HHG Amplitudes: 
Many-molecule vs. One-molecule Signal 

It is interesting also to consider the total amplitude 
Atot.(ty °f HHG emission from all the molecules inter- 
acting with the (probe laser) field. This is given by the 
coherent sum of the individual amplitudes emitted by the 
molecules at the positions Xj for all j — 1, 2, 3, ■ • ■ A/", 
where Af is the number of molecules in the interaction 
volume, or 

M 



A tot .(Q) = J^A,(n) 



e(e 



-liKn-nk^-Xj 



x \ -2iri - nw)T iti (n) I (27) 

I n ) 2 

where we may identify the basic HHG transition ma- 
trix element for the emission of the nth harmonic per 
molecule as, 



(28) 



It can be seen from Eq. (|27f that the nth harmonic emis- 
sion amplitude in fact factorizes into two parts, the first 
factor corresponds to the sum of the macroscopic space 
dependent phases associated with the random positions 



Xj of the C.M.s of the molecules in the interaction vol- 
ume, and the second factor corresponds to the funda- 
mental "one molecule" emission amplitude, independent 
of the position of the C.M.s of the molecules. The macro- 
scopic phase factor is explicitly given by 



{•••}> H E 



-j — iiKa—nk^-Xj 



(29) 



For a large number of molecules in the interaction vol- 
ume, Af ^> 1, the phase factor oscillates greatly and thus 
tends to average out to zero, except when the condition, 



(Kr 



nku) = 



(30) 



is fulfilled; in that case it yields the phase sum = Af. 
It is readily understood that the condition ([3U|) corre- 
sponds exactly to the momentum conservation between 
the final momentum of the emitted harmonic photon, 
hKii, and the sum of the momenta of n laser photons, 
nhku- This is the phase-matching condition in the for- 
ward direction (37j . 

The probability of emission of the harmonics is given 
as usual by the absolute square of the total amplitude 
Eq. (|2"T|) . Under the phase-matching condition the latter 
is therefore coherently amplified by a (generally large) 
factor of Af 2 . This is also the origin of the quadratic 
pressure dependence of the high harmonic signals, as well 
as their unusual strengths, that had been found in the 
very first experimental observations (e.g. (38l. l39j). 

The second factor {...}2 gives the fundamental "one- 
molecule" quantum emission amplitude. We also note 
that if the absolute probability of the harmonic emission 

is needed then the proportionality factor Nq = 



becomes essential, and that for a given polarization di- 
rection of the emitted photon, en, the projection of the 
dipole expectation value must be taken in that direction. 



D. Continuous Medium and the Phase-matching 
Function 



If one assumes that the gas molecules are distributed 
effectively continuously with a distribution function 
Afp(R)d 3 R, where p(R) is the so-called "density func- 
tion per molecule" , then one may replace the sum over 
j in Eq. (|2"9"|) by the integration over the interaction vol- 
ume. Clearly, in this case the square of the macroscopic 



phase factor, | [. 



|i| , takes the form 

.]i| 2 ^A^TiKn-nK) 



where, 

T (K n - nk u ) = 



J d 3 Xp{X)e^ Kn - nk ^ x 



(31) 



(32) 



which is the so-called phase-matching function. It peaks 
for its argument near zero (near the forward direction), 
but falls off rapidly away from it. 
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E. Coherent Elastic vs. Incoherent Inelastic 
Transitions 



and sum over the emitted photon modes (with 
L 3 J dK n J dKnKl) and get: 



Eqs. (|3"U|) and §Fj\ show, respectively, that both the 
phase matching condition (momentum conservation) and 
the frequency matching condition (energy conservation) 
in the process ought to be fulfilled simultaneously in or- 
der that the macroscopic signal to be coherently ampli- 
fied in space and time. As already noted earlier, the 
former condition leads to the directional coherence (for- 
ward propagation) of the HHG emission, while the latter 
implies the elastic nature of the accompanying molec- 
ular transitions for which the final ( "recombination" ) 
state / of the molecular system is the same as the ini- 
tial state i, with Ei — Ef. In contrast, for an inelastic 
transition, i — ► /, when Ei ^ Ef, there would be in 
general only "hyper-Raman" emissions, with frequencies 
Clif — (ji'lo — \Ei — Ef\), that are generally incommen- 
surate with the incident laser frequency or its multiple, 
or the HHG frequency, £1 = ntu. Thus the non- vanishing 
relative phase difference A</> = (\Ei — Ef | — nu) t would 
fail to stimulate the hyper-Raman transitions by the in- 
cident field, unlike the stimulated spontaneous nature of 
the associated HHG. Also the non-vanishing momentum 
difference between hyper-Raman radiation and the the 
multiple of the laser photons Kn t . — nk u ^ makes the 
former macroscopically and directionally incoherent. 

Finally, we note that the probability of the electron- 
ically inelastic processes associated with the transitions 
into the continuum (e.g. ionization) or between continua 
(e.g. inverse Bremsstrahlung [i(|) that are commensu- 
rate with the emission of the nth harmonic at the "one- 
molecule" level, will be incoherent spatially, and therefore 
would enhance only proportional to the total number of 
molecules, A/", in the interaction volume. This is in stark 
contrast to the coherent amplification of the HHG emis- 
sion at the nth harmonic, that is proportional to Af 2 . 



Differential Rate of Coherent High Harmonic 
Generation 



To derive the explicit expression for the probability of 
HHG per unit time i.e. the rate of generation of coherent 
high harmonics, we take the absolute square of the total 
HHG amplitude Eq. P7|) and divide by the long obser- 
vation time T, use a useful representation of the square 
of the delta-function (@, p. R12]), 



T 

5 2 (n-ncj)= lim — 5(n-nu), (33) 



= T lim £ 



K 



T 



at 2 f dk ^ T - nk ^ 

n J 

x / dW(nuj) (34) 



where, 

dW(nuj) = 2tt8(Q. - nuj)L 3 

x N n e n ■ D hl (ncj)\ 2 K^dKn (35) 

is the differential rate of HHG per molecule. Noting that 
the main contribution arises from the phase matching 
condition along the forward direction, we may carry out 
the mode-integrations to get: 



W(nuj) = 2tt 



T, 



(n) 



(mo) 2 



(36) 



where we have used, Kq = — , Kq = KqKq, and the 
fundamental transition matrix clement for the emission 

(n) 

of the nth harmonic, T i \ , is given in terms of the F.T. 
of the dipole expectation value D^i(nui) by: 



(37) 



for, L 3 \N n \ = 2^(nw). 



IV. EVALUATION OF "ONE MOLECULE" 
AMPLITUDE 



HHG 



Clearly the dynamical properties of the HHG signal 
are given by the rate of HHG emission per molecule, 
Eq. (f3"6")l , while the total signal is the same to within 
a proportionality constant given by the square of the 
number of molecules in the interaction volume, Af 2 , and 
the phase-matching constant f dKnJ- (Kq — nk^) that 
peaks in the forward direction. We therefore proceed to 
evaluate the dynamical signal per molecule (in a relative 
scale) as follows: (i) solve the Schrodinger equation for 
the nuclear and the electronic motions of the interacting 
laser-molecule system, (ii) construct a complete set of or- 
thonormal reference states, \i) = \Xi(t)), of the molecule, 
created by the pump pulse, (iii) determine their statisti- 
cal weights according to the one-to-one correspondence 
with the thermally occupied rotational eigenstates of the 
ensemble, (iv) calculate the "one molecule" probability 
amplitude for HHG for each member of the ensemble of 
linearly independent reference states \i), using Eq. (|36p . 
and finally, (v) obtain the (scaled) signal "per molecule" 
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by thermally averaging the probabilities of HHG emis- 
sion from each member of the ensemble of the reference 
states, using the distribution of their statistical weights. 

In the Born-Oppenheimer approximation and non- 
overlapping pump and probe pulse condition, we may 
consider the evolution of the wavefunctions of the nu- 
clear and the electronic parts separately and combine 
them together to obtain the wavefunction of the inter- 
acting system to evaluate the transition matrix elements 
of interest. 



A. Pump Pulse Interaction and Rotational 
Wavepackets as Reference States 

The nuclear rotational motion under the action of the 
pump pulse is determined by the Schrodinger equation 
governed by the partial Hamiltonian 

H N (t) + V N - L1 (t), (38) 

i.e. 

i^ JM (t) = (H N + V N - Ll (t)) $ JM (t). (39) 
at 

We first construct the fundamental set of linearly inde- 
pendent solutions of Eq. (|3"9"j) . each evolving indepen- 
dently from each of the occupied rotational eigenstates 
{| JoA/o)}- We expand it on the basis of the eigenstates 
{\JM)}, as 

* W*) = E C JM M0) (*) \J M ) e~ iEjut (40) 

JM 

The coefficients Cj J M M °\t) satisfy the system of coupled 
linear differential equations 

44m M0) (*)= E (JM\V N . Ll (t)\J'M')C^(t) 

J'M' 

r (41) 

This set of equations can be easily obtained (e.g. [41() by 
projecting on a given eigenstate from the left. In practice 
we obtain the set of the fundamental solutions \<&j m ) by 
numerical integration using the well-known Runge-Kutta 
method (34[ , starting with the following independent ini- 
tial conditions: 

CjM Mtl) (**) = f>J,Ja&M,Ma ( 42 ) 

We may note explicitly here that (a) each independent 
wavepacket-state \<&j M (t)) evolves in one-to-one corre- 
spondence with the initially occupied rotational eigen- 
state \JqMq). Taken together they form a complete set 
of orthonormal rotational wavepacket-states (linear su- 
perposition of rotational eigenstates): 

E l*W*)> <*W*)I = 1 (43) 

Jo M 



In general a gas jet of molecules in a pump-probe experi- 
ment at a finite temperature T, is not in a pure quantum 
state but rather is in a state of thermal mixture of of the 
rotational eigenstates, {| JoMo}}. We therefore introduce 
the quantum statistical mechanical device of a hypothet- 
ical ensemble of mutually independent and identical ref- 
erence molecules, each of which occupies the electronic 
ground state and the rotational eigenstates {| Jo-Mo)}, 
the latter with statistical weights p(JqMq), given by the 
Boltzmann distribution: 

p (e, J M ) = (l) e x Z P e- E ^' kT , (44) 

where 

Z P = E (2 Jo + 1) e- Ej ° /kT (45) 

Jo 

is the rotational partition function; -Ej .a/ = 
Jo (Jo + 1) hBc, for all Mo; B stands for the rotational 
constant. We shall assume for the present purpose that 
the pump pulse is not too strong so that the change in 
the occupation probability of the ground electronic state 
after the pump pulse interaction is negligible and hence 
the ground electronic state at a time t before the inter- 
action with the probe pulse evolves simply to \4> e (t)) = 
e -iE r t 1 e (0)} , where E e is the ground state energy. (We 
may assume that electronically only the ground electronic 
state |0 e (O)) is occupied initially). Thus, the linearly in- 
dependent reference states of the molecule, after the in- 
teraction of the pump pulse and immediately before the 
interaction with the probe pulse, can be written as the di- 
rect product of the nuclear rotational wavepacket states 
and the electronic ground state: 

\ Xi (t)) = |*Wi)> |0 e (i)) * = (e, JqM ) . (46) 

The reference density matrix describing the molecular 
ensemble prepared by the pump pulse takes the form: 

Pmol (e,J M ) = J2\Xi(t))p(e,J M )( X i(t)\ 

i 

= \Mt)) \$j M (t))p(j M ) 

x <*W*)I fae(*)| (47) 

where, i = (e, JqMq). The above ensemble of molecular 
states describes the effective "initial" condition of the sys- 
tem after the pump pulse, when the probe pulse arrives 
at the molecular at Xj. To avoid any possible confusion 
regarding the presence of the "mixed-state" of the en- 
semble, and the "rotational coherence" , we may already 
point out explicitly that while the ensemble is charac- 
terized by the statistical occurrence of the orthonormal 
reference states {\Xi=e,J M (*))}) each one of these states 
carries the information of the rotational coherence in- 
duced by the pump pulse, as coded in the individual 
rotational wavepackets {\®j M {t))}- Thus, when the 
thermal average of the HHG emission signal must be 
taken with respect to the probability of emission from 
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each member of the ensemble (as required by quantum 
statistical mechanics), it can not, and will not, wash out 
the rotational coherence that is present within each of 
them individually. 



B. Interaction with Probe-Pulse and Evolution of 
the Electronic State 

To proceed further, we next consider the evolution of 
the electronic state, governed by the partial Hamiltonian 



H e + V e - L2 (t - t d ) . 



(48) 



It is obtained conveniently from the knowledge of the 
electronic Green's function G e (t,t') [2(| associated with 
the above Hamiltonian and defined by the inhomoge- 
neous equation: 

- {H e + V e -L 2 (t - fd))} G e (t, t') = S(t - t')l. 

(49) 

A solution of the above equation can be written as (in 
the strong-field KFR-approximation) : 



G e (M') = -iO(t-t')J2\^ i 
x (p (*' - t d ) 



(+)\ e -iB+t 



P(t~ U) 



iE+t' 



4+> 



where j runs over all the ionic electronic states 



(50) 



with eigenvalues Ej~ , of the molecular ion and p is the 
free momentum of the electron; p(t) stands for the instan- 
taneous momentum in the presence of the field, defined 
as p{t) = (p + . The validity of Eq. ([50]) (within 

the Born-Oppenheimer and KFR approximation) can be 
verified by substituting it in Eq. (|49[) and using the com- 
pleteness relation 



(51) 



of the of Volkov wavefunctions defined by: 

(r \(t> p {t)) = e^M^^Cf (p 2 («))d«) (52) 
as well as the completeness relation of the ionic states 



El 



/,(+) 



/,(+) 



= i 



(53) 



We should note that the ionic states are generally much 
more tightly bound than the active electron in the highest 
occupied molecular orbital (HOMO). Thus in deriving G e 
above, we have further neglected the change in the ionic 
states due the interaction with the probe pulse, which we 
may refer to as "bare-ion" approximation. 



Finally, using Eqs. (|43|) and (|50|) we obtain the total 
Green's function Go(t,t') of the interacting system: 

G (M') = -i6{t-t') \^ } )\Mt-td)) 

jpJM 



x(cf, p (t'-t d )\^ +) 



(54) 



The above Green's function of the system (|54|) therefore 
holds under (a) the adiabatic Born-Oppenheimer, (b) the 
strong-field KFR and (c) the "bare-ion" approximations. 



C. The Total Wavefunction in Strong-Field 
Molecular KFR-approximation 

Combining the Eqs. (JTTJ) and (14^)) . we obtain the total 
intense-field molecular wavefunction at the lowest order 
strong field KFR-approximation: 



\Mt)) = \Xi(t)) + / dt' 

J — CO 

xG (t,t')Ve-L 2 {t' -td)\Xi(t')) (55) 

where, 

\ Xl (t)) = \Mt))\^J M (t)) (56) 
is a member of the ensemble of reference states of interest. 

D. Evaluation of the Dipole Expectation Value 

In the above we have obtained the necessary ingredi- 
ents for evaluating the expectation of the dipole operator 
Eq. (23]) explicitly. Substituting Eqs. (54]) and (55]) in 
Eq. (H, we get: 

Di,i(t) = ( Xi (t) d e y\ 1} (t))+c.c. 



dt' (<t> e (t)\ (®JoM (t)\d e 

) 

jpJM 

x |$. 7M (t))e-^ +(t - ($.7M(i') 

x(4> p {t' -t d )\(j> { + ] 

X V e - L2 (t' -t d )\$j oMo (t')) \(j) e {t')) 



(57) 



To simplify the expression further we (i) change the 
variable t' — * t' + t d , and similarly, t — > t + t d ; (ii) 
note that the free evolution of the rotational wavepacket 
after the interaction with the pump pulse at (t + t d ) 
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is $j m (t + t d ) 



$j m (td), and similarly, at 



(t' + td), ®j oMo (t' + t d ) = e- iHNt '$j oMo {td), (iii) the 
time dependence of the unperturbed initial electronic 
state at e (t + td) — \4> e ) e~ %Ei ( t+td , and similarly for 

the ionic states, <^ +) (t + t d ) = (j)^^ e~ iE t { - t ~ td \ (v) 

introduce the overlaps ("Dyson-orbitals"), 

#) = ( 4>\ +) (1, 2, ..JV e - 1)| 2, ...JV - 1, TV)) 

(58) 

and retain only the (dominant) contribution from the 
lowest lying ionic state (j — 0), to obtain: 



A,i(t) 



"* \®JoM (t) 
jJM.p 



,(0) 



X — *<*)) \<f>f) |*JAf(«) 

J — OO 

xe" l/ ''-^ (p( " )2/2)dn 
x($Jm(*')|(<#| (*>(*' -*<i) 



X 
+C.C. 



7 e _L 2 (f - t d ) |^ 0) }|$J Mo(*') 



(59) 



Or, 



= ($j M (t)\D e (t)\$ JoMo (t)), (60) 

where clearly the electronic part of the expectation value 
is given by the integral 



D e (t) 



p 

td+t 



d. 



P(t) 



' dt'e~ i ^' ( - p{u)2/2+EB)du 



x (p(t')\-F(t') • 4|4 0) )} + cc. (61) 

where, -F(t) is the probe field. Finally, by using the rate 
of emission of the nth harmonic as given by Eq. (|3"6")) . we 
obtain the dynamic HHG signal, for a pump-probe delay 
time td, 

SW(t d ,a) = 2n p{Jo)\($J M (td)\T( n Hd,<t>;a) 



Jo M 



X $ 



Jo M 



(td] 



{nujf 



(62) 



where (9, <p) = Rn is the orientation of the molecular axis 
m space and Tt\0,(f>; a) is a HHG transition operator. 

E. Derivation of the HHG operator T (n) (6>, 0; a) 

We shall now proceed to derive an explicit expression 
of the HHG transition operator Te (6,<p;a) for an ar- 



bitrary relative angle a between the linear polarization 
directions of the pump and the probe pulse. To this end, 
we first consider the most common experimental geome- 
try in which the pump and probe polarizations are chosen 
to be parallel. 



1. HHG operator: Parallel polarization a = 

We recall that for a linearly polarized probe pulse 
F 2 (t) = en F(t) cos ujt, the corresponding vector poten- 
tial is 



A(t) 



en 



cF(t) 



sin tot 



(63) 



It is convenient in this case to take the space fixed polar 
axis (z-axis) along the common direction of the polar- 
izations ei || €2 || z. To evaluate the tripe-integral over 
the intermediate momenta p in Eq. (|6ip we employ the 
stationary phase method [42j , with the stationary values 



Pst(t,t') 



1 



t-t< 



A(t")dt", 



(64) 



for which the derivative of the action S(t, t') with respect 
to t' = t — t is equal to zero. The corresponding station- 
ary value of the action is, 



S,t{t,t') = 



+- Eb j dt", 
(65) 

where, p(t) = p s t(t,t') - \A{t) and p(i') = p st (t,t') - 
-A(1/). Thus, projecting the resulting value of D e (t) on 
to the polarization direction en of the emitted harmonic 
we get, 



D e (i) 



Jo \e + i(t-t')/2 / 



s(0) 



en.r\p(t))e 



-iS st (t,t') 



p(t')\F{t').r\cl>W)}+c.c. (66) 



We may note that the first matrix element in this ex- 
pression (reading from the right to the left) corresponds 
to the "ionization" transition at time t', di on (t') = 

(p (t') \F(t').r\ 0^), whereas the last matrix element 

corresponds to a "recombination" transition of the elec- 
tron back to the same initial state at a time t, d rec (t) = 

(^e" 1 \eQ.r\p(t)^. The interval (t—f) corresponds to the 

intermediate time that the electron spends in the contin- 
uum Volkov states, between the absorption of n photons 
in the first step and the emission of the harmonic fre- 
quency f2 = mo in the last step. We have assumed that 
there was no significant depletion of the ground state 
population during the process. However, if needed, a 
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weak depletion due to ionization could be accounted for 
without difficulty by introducing in the above expression 



the exponential decay factor: 
total ionization rate. 



-(7/2)(t+t') 



where 7 is the 



2. The "ionization" and "recombination" matrix elements 

To evaluate the matrix elements of "ionization" and 
"recombination" in Eq. (|66[) . we assume that the wave- 
function of the active electron may be given by the 
highest occupied molecular orbital or HOMO. This can 
be written either in the multi-center LCAO-MO form, 
or, by transforming it into an equivalent single-center- 
MO (by appropriate transformations, e.g. (43|, 0]). It 
is useful also to note that in the latter form, it often 
suffices for problems of interaction of molecules with 
long-wavelength laser fields, and in conjunction with the 
"length gauge" , to retain only the asymptotic limit of the 
orbitals at distances away from the molecular center. 

In general, let the unperturbed MO of the active elec- 
tron of a linear molecule, in the Born-Oppcnhcimcr ap- 
proximation, be given in the body fixed frame by the 
single-center expansion: 



Ri(r)Y lm (r) 



(67) 



where C, are the expansion coefficients (at a given in- 
ternuclear separation R), Ri(r) are the radial waves of 
angular momentum I, Yi m (r) are the spherical harmon- 
ics, and to is the projection of the angular momentum 
of the active electron along the molecular axis, that is 
conserved and hence also characterizes the MO. 

Next, we transform the molecular orbitals Eq. (|67|) 
from the body fixed frame to the space fixed frame by 
using the Wigner transformation D, 



^°>(r) = 



1 



(68) 



Above D'„ 



>ro (0)e~ im * is the Wigner 
rotation matrix where (<j), 9, \) are the Euler's angles 
which define the orientation of the molecular axis to the 
space fixed coordinate frame (45[. The middle term of 
the Wigner matrix, d^ m (0), have been tabulated e.g. in 



reference [45|, |46j . The matrix element of the dipole along 
the direction of the probe pulse, appearing in Eq. ([66]) . 
then reads: 

dion ( t ') = F(t')j2cl m) Y, D U(<t>> e >x) 

I M 

x(e i P i - r \e u .r\R l (r)Y ll ,(r)) (69) 



Further, we expand e lPt ' r in spherical harmonics, 



a tp t ,.r _ 



2tt 



I'm' 

xYvmf (Pf) Yvm,' (r) (70) 



and note that in this system of axes we have, 




(71) 



Note also that the instantaneous momentum p t < can be 
either parallel or anti-parallel with respect to the direc- 
tion of the field so that 6 Pt = 0, it and <f> Pt = 0. Therefore, 
the spherical harmonics with the argument pt' can be 

simplified to Y Vm ,(p v ) = (a) 1 ' <J^±5 m > t0 , with a = 1 
for 9 Pt = and a = — 1 for 9 Pt = tt. Substituting Eqs. 
(ffil)) and OH) in E q- & we obtain (with V = k), 



d ion {t') = F(t')Y / C^ ) D l ^ Q (cl > ,9,x)/3 l on(k,m;t')(72) 



(3io n (k,m;t') = J^iiv) 1 VVk + 1) 

x (k0 \10\lfj,) I ia (t') 



(73) 



where, we have defined the radial integrals (to fixed) by 



Iia(t')= J^Apt^RTMr-^r^dr (74) 



Using the Slater orbitals representation of the single cen- 
ter radial functions, Ri{r) — r n ~ 1 e~ ar , the radial inte- 
grals (I's) appearing in d ion (E g (1721) ) can be evaluated 
explicitly by using the formula [47| 



e-^J^i^x^dx = 



sj{p? + /? 2 )"+T(l/ + 1) 
+ 1 — fl + V 



xF 



v+1; 



/? 2 



a 2 + (i 2 



(75) 



where F (a, 6, c; x) is a hypergeometric function. Note 
that since the argument x = a2+l3 2 < 1 the hypergeo- 
metric function is guaranteed to converge for all values 
of a, b, and c. For ionization step, the radial integration 
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reads 



V(PB+P 2 t>) h+Zc/PB+4 T(h + I) 
k + Z c /p B + 4 k - Z c /p B 



xF 



k + 



pI 



2 Pi + Pi 



1,2 T{l l + Z c /p B + 2) 

V(p 2 B +pi) M/pB+2 nh + h) 

k + Z c /p B + 2 k - Z c /p B - 2 



xF 



k + 



2'PB+Pf 



(76) 



Exactly the same expressions hold for the radial integrals 
appearing in the recombination dipole d rec (l r ,m;t) (ex- 
cept that U is changed to l r , and t' to t, in Eq. {76])), 
throughout. 

The Clebsch-Gordan coefficient in Eq. ([73]) implies 
that only the /x = term and I = li ± 1 terms survive in 
the sums and we get for the ionization dipole, 



with 



(77) 



2tt 



V(2X 



x ((ia)^ 1 (k + l)I iai+1 (t') 

+ {iaf- l l i I UM _ 1 {t')) (78) 

In the above we have introduced the angular momentum 
notation U for the initial bound state in the "ionization" 
matrix element (and l r , in the "recombination" matrix 
element). 

We may assume that the emitted harmonic is observed 
with its polarization along the same direction as the 
probe pulse polarization. (There is no difficulty, except 
lengthier algebra, to obtain the expression for the polar- 
ization direction orthogonal to it, but the former would 
give the dominant contribution under phase matching 
condition. Following an analogous calculation as above 
we get the "recombination" matrix element as: 

d rec (t) = ]T Cf^* Dfofa e, X )Prec (lr, mj f ) (79) 



with 



p rec (l r ,m;t) 



2tt 



x ((ia) l - +1 (l r + 1) Il r ,l r+ l{t) 

+ {ia) l ^ 1 l r I lr , l ^ 1 {t)) (80) 



Substituting Eqs. ( I77ll78|) and ([791180]) in Eq. ([66]), we 

obtain 



D e (t) = i£4«(0)4n(0) 

li,l r 

x \($ n) *c£ l) M e (t)+c.c 
where we have used the relation 

D% (4>, 6, x)D l ^ (<t>, o, x) = 4 m {6)d l t m (9) 
and defined the radial integral, 



(82) 



rt d +t , N 3/2 

MJi) = i dt\ t -, -r-r. 



x / 3 rec (Z r ,m;0e-^ t(t ' t ' ) 
xF(t')f3 ion (k,m;t') 



(83) 



Next, by integrating over t' , taking the Fourier trans- 
form with respect to t, we obtain (cf. Eg. ([28]) ), the HHG 



operator T t 



(«), 



; 0) for the nth harmonic generation: 



Tj")(0,0;O) = y/2n(nu)D e {ntj) 

li,lr 

Xa'" ] (l r ,k;m) 



(84) 



where, (l r ,li;m) is given by the nth Fourier coeffi- 
cient of D e (t): 

4™ } (lr,k;rn) = 

F.T.[cl™ ) *c\fM e {t) + c.c'\ (nu) (85) 

we obtain 



Next, by substituting Eq. ([HI]) in Eq 
the rotational matrix elements, 



$J M (t d ) 4n(0)4n(0) *J M (t d )), (86) 



which can be evaluated directly by using the tabulated 
values of the d l Qm (8) given by elementary trigonometric 
functions (see, Tab. [I]). Alternatively, we may first com- 
bine the product 



47T 



-Y, 



li,m\ 



and we obtain 



2k + 1 

53 (-l) m ^,^,-m,m|L,0) 

L=|i r -Ii 

x (Z r ,Zi, 0,0|L,0)P L (cos 6) (87) 



rj">(0,<£;O)= ^ a'i" J (Z r J i ,L;m)P i (cos( 

h,l r ,L 
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where, 

a'^J (l r , h, L; m) = y / 2ir(nuj)a z 7 l > (l r ,k;m) 

x(-l) m (l r ,k,-m,m\L,0) 

X (l r) l h 0,0|L,0) . (89) 

Thus, the expectation value of the transition operator 
with respect to the rotational wavepacket can be obtained 
more elegantly in terms of the Legendre polynomials mo- 
ments: 

(Pl)j oMo (*«0 = (*J M (t d )\ P L (9) |$ JoMo (U)) (90) 



Finally, by substituting the above relations (Eq. ([55]) ) 
in Eq. (|62j) . and taking the statistical average over the 
ensemble of the emission probabilities from the ensemble 
of rotational wavepackets, we obtain the HHG signal (i.e. 
the rate per unit time of generation of the nth harmonic 
per molecule) in the special case of parallel polarizations 
(cf. 0): 



SW(t d ;0) = 2tt £ p(J ) 

Jo M 



a 'zz Q-rih, L;m) 



X (PL)j oMo (t d ) 



(91) 




Figure 3: A schematic diagram defining: molecular axis R, 
electron position r, pump polarization ei, and probe polar- 
ization 62; a is the operational laboratory angle. 



F. General Polarization Geometry: Arbitrary a 

So far we have assumed that the pump and the probe 
polarizations are parallel and that they point along the 
space fixed polar axis z. In the general case, we may 
define, without loss of generality, the relative angle be- 
tween the polarizations, a, to lie in the (z — z' — :r)-plane 
(cf. Fig. [3]). From the figure, it can be seen that we 
simply need to re express the direction of the molecular 
axis, (0,0), given with respect to the pump polarization 
ei || z), in terms of the direction (#',</>') with respect to 
the probe polarization direction €2 \\ z. This is readily 



achieved by simply replacing cos 9 — > cos#', and using 
the well known relation 

cos 9' = cos a cos 9 + sin 9 sin a cos <fi (92) 

or the vector addition coefficients and the addition theo- 
rem 

Pl(cos9') = ^—^Y LM (9, 4>)Yl, M {a, 0). (93) 

Thus, we obtain the general expression of the HHG op- 
erator for any a [20| : 



k,lr 

xaW (l r ,k;m) . (94) 



Or, 



^ a' z " (l r , k,L;m)P L (cos9') 

l T ,U,L,M 
LM l r ,h 

x^jY L , M (d,cj>)Y^ M (a,0) (95) 



where a' zz (l r , k, L; m) is given by Eq. (f89|) . 

It is should be noted that, in general, if the molecular 
orbital coefficients were assumed to be complex, 



C, 



(m) 



(96) 



then we should rewrite the dynamic parameters 

Ct'z) {lr, h; m) 
(Eq.©) as: 

(l r ,k;m) = F.T. 



cl m> ci m) M e (t)+c.c. 



(m) 



a 



cos cpuir-u) \. 



sm( Phl r v l r ,U, m 



(97) 



where, 



4> h i r = (ct> h -<p tr ). (98) 
M e (t) is given by Eq. (f83|) . and we have defined, 

= F.T. [Re {M e (t}}] (ncu) (99) 



-(") 



r,(«) 

l r ,li,m 



V 



F.T. [Im{M e (t)}] (nw) 



(100) 



Thus, finally, we can express the general transition ma- 
trix element for the nth order harmonic as an expansion 
in Legendre polynomials in cos a, and the corresponding 
Legendre moments of the time-dependent axis distribu- 
tion of the molecule: 



14 



cos 



(4>U,i r ) u 



li,l r ,m 



2 | C M|| C M| 



TW (t d) a) - ($./ Afo fa) I rW(S, 0; a) |$ JoMo 
= ^ \/ 27r(na;) 

-sin (<j>t t i r ) vfcj 

x(-l) m {l r ,k,-m,m;L,0) 
x (Z r , Z i5 0,0;i,0> 

x (PL}j oMo (t d )P L (cosa) (101) 

where, we have taken the expectation value of the 
HHG operator with respect to the rotational wavepacket 
|<I> (td)), to obtain 



4~ 



2L + 1 



(*JoM (^)l>LM(e, 
= <^>,/„M (td) 



>>)\$J M (td)) 

P L (cosa)6 M ,o 



(102) 



This follows from the observation that the magnetic 
quantum numbers of all the rotational eigenstates in the 
individual wavepackets have the same value M . It holds 
when the interaction operator V e -L 2 (t) does not depend 
on the azimuth angle of the molecular axis in the body 
fixed frame. In the above expression, we have also used 

the relation 



2L + 1 



L0[ 



) = Pl (cos 9), and an analo- 
gous relation with respect to the angle a, to simplify. 

It is useful to note also that, if the orbital expansion 
coefficients are real, as is often the case, then 4>i it i r = 
(0, 7t) and therefore the quantity in the square brackets 

in Eq. (| 1 1 [) simplifies to [...] = cos (</f; ij ; r )u| ra ! only. 



Table I: Explicit form of d' 0m (6) required for evaluating Eq. 

dm) mm 
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N 2 (m = 0) 




o 2 


(m = l) 





1 








2 


| (3 cos 2 6»- 1) 




[z 

V 2 


sin 8 cos 6 


4^(3 


- 30 cos 2 + 35 cos 4 s) 


4 


sin 9 


cos 0(3- cos 2 6») 



APPLICATIONS TO DIATOMIC 
MOLECULES N 2 AND 2 



A. Parallel Geometry, a — 0: Elementary 
Expression of Tj n) (0,a = 0) for N 2 



N2 has (j g symmetry, and we approximate its MO 
by the asymptotic approximation from single center 
molecule (Eq. JB7J) with m = and I = 0, 2, 4 H 
whose angular coefficient are given in table [TTJ The radial 
part of electronic wave function is given by 



Mr) 



(104) 



G. A General Formula for the HHG Signal 



Thus, finally, we substitute Eq. (flUTj) in Eq. (JHU) and 
obtain the desired general expression (cf. [20j) for the nth 
harmonic signal from a linear molecule, for any value of 
td and a: 



S {n \t d ,a) 



^2 PJoM„2tt 

Jo M Q 

■ J. -M 

- sm ipt f t r vi \. 



COS (pi I u 



l r ,li ,m 



L.l r ,L 



2 | C M||(7M| 



J i T w^i 
x(-l) m (l r ,l i: -to, m;L,0) 
x (l r ,U,0,0;L,0) 

x ( p l).JoMo (td) Pl (cos a) 
(niu) 3 



(103) 



We may conclude the section by noting that for the 
special case of parallel polarizations, a — 0, Pl(cosO) = 
= 1, Eq. (|103[) correctly goes over to the signal 
obtained for that special case, Eq. (p?Tj) (cf. (l9j). 



with 77 = Z c /pb'i Z c is the core charge and — y/2 \Eg\ 
with .Eg is binding energy. 



Evaluating Eq. (|84|) for to = and k,l r — 0, 2, 4 give 
us the HHG operator for N2 



Zi,Z r =0,2,4 

(105) 

Using the expressions for the reduced rotation matrices 
from Tab. Uand simplifying, we may rewrite the operator 
as a sum of powers of cos 2 9 only, 



T^(9) = ^2M^7) [4" ) +6 1 " ) cos 2 + 6 2 " ) cos 4 ^ 

+b^ n) cos 6 6 + b { 4 n) cos 8 9} (106) 



15 



where the coefficients b^s reads, 



ai™;(0,0;0)-ia^(2,2;0) + ^(" ) (4,4;0) 
-i(4" ) (0,2;0)+^)(2,0;0)) 
+|(aW(0 ) 4;0) + oW(4,0;0)) 



k(«) 



--(oW(2,4;0) + oW(4,2;0) 
-|sW(2,2;0) 

+|(aW(0 ) 2;0) + oW(2,0;0) 
-H(oW(0,4;0) + oW(4,0;0) 
-^(a^(2,4;0) + 4™)(4,2;0) 
^(aW(0,4;0) + aW(4,0;0)) 
(aW(2,4;0)+aW(4,2;0) 



125 

~ 16~ 
105 

16~ 

(4,4;0) 



(oW(2,4;0)+aW(4,2;0) 



(107) 



Thus, by using Eq. (p2^1 . the nth harmonic signal for N 2 
becomes, 



3 3 



S^(t d ) = C^^4;»((cos 2 ^)(cos^ 

j=0 j'>j 

= c{4 l) +4 l) ((cos^))(t d ) 

+<£> ((cos 2 Of) (t d ) + c$ «cos 4 #)) (t d ) 
+ --- + ci ? P((cos 8 0) 2 )(t d )} (108) 

where C = (y/2wnw) 2 2tt^I = (2tt) 2 (if) 3 . The coef- 
ficients cj™? are related to as follows 



„(«) 



2/ i V:(^" ) ^ l) * 



for j = j' 
forj 



(109) 



The leading two terms of the signal for N 2 , Eq. (|108p . 

(n) 

consist of a constant term proportional to Cq , that arises 
from the leading angular momentum term / = of the 
active molecular orbital of N 2 , and a term proportional 
to the second moment ((cos 2 #)) (td) that corresponds to 
the usual "degree of alignment" A (td). We may note in 
passing that the above result does not support a recent 
model calculation [53|, [54[ that emphasizes that the lead- 
ing contribution for HHG signal from N 2 arises from the 
fourth moment (cos 4 (9^; that would require, for exam- 
ple, dropping the basic contribution of the I = term i.e. 



Table II: The molecular properties used in this work. I p 
is adiabatic ionization potential, B is rotational constant of 
molecule, Qy and a± are parallel and perpendicular polar- 

izability, and C^'s are angular coefficient of the electronic 
wave function. 





N 2 


o 2 


Ref. 


HOMO 


a g , m = 


TY g , m — 1 


[22, 50] 


Ip (eV) 


15.58 


12.03 


[48] 


B (cm" 1 ) 


2.0 


1.4377 


w 


-II (A 3 ) 


2.38 


2.35 


[52] 


ax (A 3 ) 


1.45 


1.21 


[52] 




2.02 




[48] 


Q( m ) 


0.78 


0.62 


[48] 




0.04 


0.03 


[48] 



b ( o l) in Eq. (fl06|) - for the HHG operator for N 2 - that of 
course would not be justifiable due to the a — symmetry 
of its active orbital. 



B. Parallel Geometry, a — 0: Elementary 
Expression of ri n) (0;a = 0) for 2 

2 has Tig symmetry, and thus we approximate its MO 
by the asymptotic approximation with m — 1 and I = 2, 4 
p8l . |49| whose angular coefficient are given in table [TTJ 
The HHG operator (Eq. 1(51))) for 2 reads 



T^(9) = ^(nuj) 

x E d^(9)aif(l r , 1^1)^,(9) (110) 



L,l r =2A 



By using the expressions for the reduced rotation ma- 
trices from Tab. U and simplifying, we may rewrite the 
operator as a sum of powers of sin 2 9 cos 2 " 9 only 

Tj") (9) = ^2k(uu) \b[ n) sin 2 9 cos 2 9 + b 2 n) sin 2 9 cos 4 

(111) 



+b ( 3 n) sin 2 9 cos 6 I 



where b^ J -coefficients are given by 

b^ = ^)(2,2;l) + ga(«)(4,4;l) 
3 /15 

105 



u 2 — 



(aW(2,4;l) + fl W (4,2;l)) 
S^(4,4;l) 



+|VT («i" ) (2,4;l) + ^)(4,2;l)) 



245. 
16 



5^(4,4;!) 



(112) 



16 



Finally, substituting operator expression (Eq. (j 1 1 1 [) ) in 
Eq. (f6"2")l we obtain the nth signal of O2 reads 

3 3 

5fW (t d ) = c ^ c tr (( sin2 cos2 ' ? e ) ( sin2 cos2/ ) ) 
./ 1 />./ 

= C{cff ((sin 2 cos 2 0) 2 )(t d ) 

((sin 2 6> cos 2 0) (sin 2 cos 4 0)) (t d ) 
+ -" + c 33 ) ((sin 2 0cos 6 0) 2 ^ (f d )} (113) 

Above, coefficients cj", are related to coefficients of 
Eq. pT2| through Eq. ([109]) . 

We note that, unlike in the case of N2 considered above, 
now there is no constant leading term in the signal for O2, 
Eq. (|113p . This is a consequence of the ir — symmetry 
of the active orbital for O2, which does not permit the 
lowest I = angular momentum component for its active 
orbital. 



C. Arbitrary Relative Polarization Angle a: HHG 
Signal 

We now consider the signals for N 2 and O2 in the gen- 
eral case in which the probe and the pump polarizations 
make an arbitrary angle a between them, as shown in 
Fig. [3l Unlike alignment angle 0, the pump-probe angle 
a can be controlled in the laboratory and may provide 
a possible control of HHG of molecule. To obtain the 
signals in terms of elementary trigonometric functions in 
this more general case, we refer to Fig. [3] The direc- 
tion of the molecular axis is now denoted by (0', 4>'). The 
same expression for the signal as in the parallel case now 
holds in terms of the primed angles. The HHG signal 
(Eq. (fT08|) for N 2 and Eq. (|TTg|) for 2 ) for arbitrary 
angle a now can be written as: 

S™(t d ;a) = C + cg[> ((cos 2 6")) (t d ) 
+0^ ((cos 2 0>) 2 ) (t d ) 

+ --- + 4 , 4 ) ((cos 8 0') 2 )(M} (114) 

for N2 and 

S^\t d -a) = C{4 l) ((sin 2 0'cos 2 0') 2 )(^) 

+4° ((sin 2 0' cos 2 0') (sin 2 0' cos 4 0')) (t d ) 
+ -" + c 33 ) ((sin 2 0'cos 6 0') 2 )(i d )} (115) 

for 2 . Above, «/(0')» (t d ) 

Ej oMo p(Jo){$J M (td,0)\f(0')\®J M (td,0)) is an 
expectation value of function f(0') given in probe frame 
but evaluated with respect to the rotational wave packet 
obtained in the pump frame. Before evaluating the 



above integral, it is convenient, therefore, to transform 
the HHG operators in the variables (0 1 , </>') in the angles 
(0, tf>) of the pump-frame (i.e. with the z along the 
pump polarization). This can be done by the simple 
transformations, given by Eq. (|92[) , where <f> is the angle 
between plane of molecular axis - pump pulse and plane 
of pump - probe pulses. The elementary expression 
for the expectation value of the alignment operator 
A(t d ;a) — (cos 2 0') in the case of non-zero angle a 
reads: 

A(t d ;a) = (cos 2 0') 

= ^cos 2 a - i sin 2 o \ (cos 2 0) + i sin 2 a 

+- sin 2 a ((sin 2 0e 2i ' p ) + c.c.) 
+-sin2a((sin0cos0e^) + c.c.) (116) 

where (sin cos 0e ±l ^) couples the J' states with A J = 
0, ±2 and AM = ±1 whereas (sin 2 6>e ±2 ^) couples the 
J' states with A J = 0,±2 and AM = ±2. We note 
that for the linearly polarized pump pulse of the present 
interest, the interaction Hamiltonian is proportional to 
cos 2 0, which is independent of M in the space fixed 
pump- frame. Thus the M -quantum number of the ro- 
tational wave-packet remains constant, or M = Mo, 
throughout the evolution. Hence, the expectation val- 
ues of (sinflcos^e* 1 *) and (sin 2 0e ±2% ^ vanish and we 
get, 

(cos 2 0') = i (3 cos 2 a - 1) (cos 2 0) + ~ sin 2 a (117) 

In a similar way, we obtain the expectation value for 
higher order 

(cos 4 0') = - (35 cos 4 a - 30 cos 2 a + 3) (cos 4 0) 
8 

3 

+ - (- 10 cos 4 a + 12 cos 2 a- 2) (cos 2 0) 
8 
3 

+ -sin 4 a (118) 
8 

We note in passing that for a = 0, (cos 2 0') in Eq. Q117[) 
and (cos 4 0') in Eq. (|118[) reduce to (cos 2 0) and (cos 4 0), 
respectively. 

Thermal averaging Eq. (I117[) gives us the "degree of 
alignment" or the alignment-moment: 

A(t d ,a) = - (3 cos 2 a - 1) ((cos 2 6)) (t d ) + - sin 2 a 

( 11Q ) 

which also appears in the second leading term of the sig- 
nal for N 2 , for arbitrary angle a, (see, Eq. Q114p ). Squar- 
ing and taking the thermal average of Eq. Q117P gives us 

((cos 2 6") 2 ^), which is the third term of HHG signal of 

N2. The thermal average of Eq. (|118p gives us ((cos 4 0')) 
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Figure 4: Comparison of the experimental [13 , 1551 ] and the 
theoretical dynamic 19th HHG signal for N2; pump intensity 
I = 0.8 x 10 14 W/cm 2 , probe intensity I = 1.7 x 10 14 W/cm 2 ; 
duration 40 fs, wavelength 800 nm, temperature 200 K. 



Figure 5: Comparison of the experimental [13|, [5E 
theoretical 19th harmonic dynamic signal for O2; 



and the 



pump in- 
tensity I = 0.5 X 10 14 W/cm 2 , probe intensity I = 1.2 X 10 14 
W/cm 2 ; durations 40 fs, wavelengths 800 nm, temperature 
200 K. 



which appears in the fourth term of HHG signal of N2. 
The difference of Eq. (fTT7|) and Eq. (fTTg)) gives us 

(sin 2 & cos 2 0') = - (-35 cos 4 a + 30 cos 2 a- 3) (cos 4 (9) 



(30 cos 4 a - 24 cos 2 a + 2) (cos 2 0) 



1 



H — ( -3 sin a + 4 sin a) 



(120) 



Squaring and then thermally averaging Eq. (|120[) yields 
the leading term of HHG signal of O2, given by Eq. 
(j!15p . In a similar way, we can explicitly exhibit the 
a-dependence of the higher order terms in the signal for 
2 , Eq. pT5[) , as well. 



VI. RESULTS AND DISCUSSIONS 
A. Signals in the Time Domain 

We now apply the theory to analyze the observed HHG 
signals from the diatomic molecules, N2 and O2. In typi- 
cal recent experiments (e.g. (HI H2, E3, H3] ) an ensemble 
of N2 or O2 molecules is first set into free rotation by 
a femtosecond pump pulse. The HHG signals were de- 
tected by monitoring the emission due to a second more 
intense femtosecond probe pulse, that was delayed with 
respect to the first by successively increasing the time 
intervals, td, in the picosecond domain, between them. 

In the experiments for N2, for example by Miyazaki 
et al. [H, [55j], a peak pump-intensity Ii = 0.8 x 
10 14 W/cm 2 , a peak probe-intensity I2 = 1.7 x 
10 14 W/cm 2 were used; the central wavelength A = 
800 nm and the pulse duration t =40 fs were kept the 
same for both the pulses. For the experiment with O2, 
the harmonic signal was measured in a similar fashion 
for h = 0.5 x 10 14 W/cm 2 and I 2 = 1.2 x 10 14 W/cm 2 ; 
the other parameters were kept the same as in the case 



of N2. For the purpose of a direct comparison, our cal- 
culations were performed for the same parameter values 
as in these experiments [H, |55T |. In Fig. 3] and Fig. [5[ 
we compare the calculated HHG signals as a function of 
td for N2 and O2, with the experimental data obtained 
for the 19th order harmonic. The effective ensemble tem- 
perature was taken to be T= 200 K, that was estimated 
from the matching of the peak position of the spectral 
distribution with that of the Boltzmann distribution as 
suggested first in It can be seen from Fig. [5]that the 
experimental data for N2 show the "revival" phenomenon 
with a full revival period T rev — 8.4 ps (which is consis- 
tent with the rotational constant of N2 (cf. Tab. [TTJ) ) 
as well as a i and a | fractional-revival. The observed 
signal for O2 shows, in addition to the full revival (pe- 
riod for O2 is T rev = 11.6 ps) and two fractional revivals 
similar to the two seen for N2, an additional i revival. 
The calculated signals can be seen to follow the same se- 
quence of the full and the three fractional revivals as seen 
in the experimental signal. We note that these observa- 
tions for N2 and O2 are also consistent with the data of 
Itatani et al. [ll[ and Kanai et al. 

To understand the similarities and the differences be- 
tween the signals for N 2 and O2, we use the analyti- 
cal results of the present theory below. The properties 
of the HHG signal of N 2 are governed by Eq. ()108|) . 
The first term gives a constant background. The sec- 
ond term ((cos 2 #)) (td) is the dominant dynamic term 
and makes the signal to mimic the "degree of alignment" 

A(t d ) = {(cos 2 9)) {t d ). The third term (|(cos 2 6) 2 ^ (t d ) 

give unequal maxima and minima i.e. the difference 
between the maximum signal and the average signal 
is greater than the difference between the average sig- 
nal and the minimum signal. Furthermore we point 
out that at a lower initial temperature, the valley of 

(^(cos 2 d) 2 ^ (td) that occurs, for higher temperatures, 
at the jT rev revival, can split into two valleys, due to 



this term, and thus the third term can strongly affect 
the HHG spectrum, as can be seen in the experiment 
by Itatani, et al. [Til . [lil ]. Another earlier puzzle re- 
garding its dynamic signal observed was the failure of 
the alignment measure A(td) = ((cos 2 9)) (t d ) to ac- 
count for the dynamic HHG signal for O2, observed by 
Itatani et al. [11]. In fact, Itatani et al. found that 
their data behaved more closely to the expectation value 
B(t d ) = ((sin 2 26)) (t d ). From Eq. (fTT3|) it can be seen 
that indeed the leading term of the signal for O2 is given 
by ((sin 2 0cos 2 e) 2 ) (t d ) = i ((sin 2 29) 2 ) (t d ), which is 

directly proportional to the observed signal. Moreover, 
the present theory also predicts that there ought to be 
modifications to this result due to the higher order terms 
in Eq. (| 1 13|) . In contrast to N2 there is no significant 
difference between minima and maxima for O2, since all 
terms in Eq. (|73[) have similar minima and maxima. The 
present theory also predicts that there ought to be mod- 
ifications to this result due to the higher order terms in 
Eq. (|113| . In contrast to N2, however, there should not 
be a significant asymmetry between the size of the max- 
ima and the minima for O2, since the terms in Eq. (|73|) 
have similar maxima and minima. The can be, however, 
quantitative contributions from the higher order terms 
predicted by the theory. In fact, as mentioned earlier, 
Kanai et al. [ll| found empirically that their experimen- 
tal HHG signals for N2 and O2 demanded heuristic in- 
troduction of operators involving higher orders of cos 9 
functions, or Legendre polynomials, as the dynamic sig- 
nal could not be well expressed in term of ((cos 2 9)) (t d ) 
only for N 2 , or ((sin 2 29)) (t d ) only for O2. In fact, the 
present theory provides an ab initio derivation of the de- 
sired general expansion of the HHG signal in terms of the 
moments of the Legendre polynomials Eq. (|103[) and/or 
of the powers of cos 2 9, e.g. Eqs. (fTT3l) and (|T08|) . 

A related characteristic of interest first observed by 
Miyazaki et al (e.g. [551 ]) is the appearance of extra se- 
ries and lines in the Fourier spectrum of the dynamic 
HHG signal for both N2 and O2, that are Raman forbid- 
den. These extra lines can not be attributed to A (t d ) = 
((cos 2 9)) (t d ), for N 2 , or to B (t d ) = ((sin 2 29)) (t d ) for 
O2. It will be seen below that the F.T. of the higher or- 
ders terms of Eq. (fT08|) for N 2 and of Eq. (TTT31I for 2 , 
given by the present theory can consistently account for 
their appearance. 

Comparing the expressions for the signals for N 2 and 
O2 and directly calculating the cffl coefficients in the re- 
spective signals it is found that the signal for N 2 is much 
stronger than that for O2, as also observed experimen- 
tally Hi. 

We may briefly discuss here the dependency of the 
HHG signals on the initial temperature, an example of 
which is shown in Fig. [5J It can be seen from the figure 
that the lower initial temperature gives a greater ampli- 
tude of revival. This may be understood as follows. A 
lower initial temperature gives a lower value of the max- 
imum of the statistically occupied Jo levels and hence 
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Figure 6: The dependency of HHG signal on the harmonic 
order and on the initial temperature of the gas. Solid and 
dashed lines are for the 19th and the 21th order, respectively. 

also a lower value of the maximum initial value of Mq, 
than at a higher temperature. As a result of interac- 
tion with the linearly polarized pump pulse (quantization 
axis along the polarization axis) at a given intensity, each 
wavepacket that evolves from a given initial | Jo Mo) state, 
can couple to the higher levels J' > max. (Jo) but can 
not raise the initial maximum value of Mq. Therefore, 
for a given intensity, the ratio of J' to Mq is higher for 
a lower temperature, and as a consequence the degree of 
alignment A (t d ) = ((cos 2 0)) tends to be also higher, im- 
plying that the molecule becomes more strongly aligned 
during a revival. 



B. Rotational Revivals: Periods and Phases 

If a linear molecule has a permanent dipole moment 
(e.g. hetero-nuclear diatomics), then the interaction 
Hamiltonian of the (pump) laser with the molecular 
frame depends on the first power of cos 9, where 9 is the 
angle of rotation of the molecular axis with respect to 
the laser polarization axis. In contrast, the interaction 
with the polarizability of the molecule (e.g. for homo- 
or hetero-nuclear diatomics) depends on cos 2 9. Thus in 
general the interaction may contain the operators cos 71 9 
with n = 1 and/or 2. Then in either case, the rotational 
wavepackets created by the later can be written in the 
form: 

$J M (t) = ^2 Cj 0+n j >Mo (t) 
j=0,l,2,3— 

x e-^ Ej o+^\J 0+n j,M ) (121) 

This can be obtained, for example, from a consider- 
ation of the perturbative solution of Eq. (f3"5|) in suc- 
cessive power of the interaction Hamiltonian, and noting 



that the rotational eigenstates couple either by P2(cos9) 
(in the absence of a permanent dipole moment) with a 
minimum (non-zero) n = 2 or by Pi (cos 9) and P2(cos#) 
(in the presence of a permanent dipole moment) with a 
minimum n = 1. It can be readily understood from the 
well known properties of the vector addition coefficients 
that appear in the integration over the product of three 
spherical harmonics (cf. para below) that the expecta- 
tion value of the iVth cosine-moment with respect to a 
rotational wavepacket at a time t — t d , takes the form: 



( C0S ^.JoMo^) = (*JoMo(t d ) 



COS 



N 9\<S> 



J M ( 



(*«0> 



= E E E 

s p=-aj=0,l,2,3— 

f~,J M * (± \ (~iJ M a (f \ 

* ^J +nj+p,Mo \ d > W +ni,M K L d) 

x «s {Yj Q +nj+pMa\ ^s,0 |X/ +nj,M ) 
X exp (Bj 0+n j +p - Ej Q+nJ )t d 

(122) 

where, the integers s and p have the same parity (even 
or odd) as the parity of N. This follows from the fact 
that cos^ 9 can be expressed as a linear combination: 
cos^ 9 = J2 S a sPs(cos9), for all s up to N, and since the 
matrix elements {Yj 0+n j, iMo \ Y Sl0 \Yj n+n] ,M„) = 0, unless, 
Jo + nj' = Jq+ nj + (p — s) > 0, and J + n 3 + (P — s ) + 
Jq + + s = even. Thus, the phase of each individual 
term of Eq. (|122p , for any given value of the integers 
Jo,n,N,j, is given by 



1 



(-Ejo+nj+p — Eja+n 3 ) td 

i2irhBct (±2pJ ± 2npj + p 2 ± p) 



2tt 



t d {hBc)2 



Jo + nj + 



p +p 



(123) 



where, we have used E, 



J,M 



J (J + l)hBc; B is the 
rotational constant, T rev = is the rotational period, 
and /i = 2nh. We note first that the quantity in the 
last parentheses above is an integer, independent of the 
value of j and Jo- We note that the maximum value 
of s or p above is N. The phase difference (Eq. (|123|) ) 
therefore equals to an even or odd multiple of w, or odd 
multiple of ^, depending on the parity of the groups of 
rotational states. Therefore, the shortest time period for 
which the phases of all terms or all terms within a parity 
group become equal in Eq. (|123|) . and hence coherently 
enhance the signal, is clearly 



T ■ — - T 

± min i\ t 

nJy 



(124) 



For times between the successive coherent enhancements 
or "revivals" , the individual phases in Eq. (|123p disperse 
away from one another and the revival peaks tend to 




4 6 

Delay time (ps) 

Figure 7: Revival structure of the moment ^sin 2 2#)) (td), 
for the case of O2, on the beat frequency: for all A J retained 
(panel a), with A J = only (panel b), with A J = ±2 only 
(panel c), and with A J = ±4 only (panel d). It is clear that 
the transition with A J = ±4 has the lowest fractional revival 
at -gTrev, and has the shortest period. 



be washed out by destructive interference, and the HHG 
signal reduces to the average or the back-ground level. 

We may summarize the above result as a "revival the- 
orem" : If the laser-molecule interaction Hamiltonian is 
characterized by the lowest power n of cos™ 9, withn = 
lor 2, 9 is the rotation angle, and if the highest dis- 
cernible (numerically significant) moment in the expres- 
sion of the signal is ($j m (*)| cos^ 9 \$j M {t)), N>1, 
then the experimental signal would exhibit as many as 
n x N revivals within a full period T r = 57^7 B is the 
rotational constant. Inversely, by counting the number 
of fractional revivals in the observed HHG signal, one 
may determine the highest order, N, and hence also the 
significant "cosine moments" (up to the order N) that 
would be necessary to fit the observed signal. We may 
note that the above theorem covers the well-known cases 
of fractional revivals discussed earlier (5f| [H, as 
special cases. 

For homonuclear diatomic molecules with no perma- 
nent dipole moment, the lowest order pump pulse in- 
teraction is due to the polarizability tensor with n — 
2. Thus for the standard alignment moment, A (td) = 
((cos 2 6*)) (td) with N — 2, we get the lowest fractional 
period Ti = hT r , and the subsequent two fractional re- 
vivals Ti , T| (defined analogously) and the full revival 
at T r , with in a period. Thus the presence of the high- 
est significant fourth cosine-moment with N = 4 would 
show the lowest = | revival, plus the subsequent 
six fractional revivals at (4, |, i, |, | , J) T r , within a 
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Figure 8: Weak components of the fractional revivals in the 
dynamical alignment signal for the case of O2. Contribution 
from ((cos 6 #)) {td), near T rev /12 (upper panel), and from 
((cos 8 #)) (t d ), near T rev /16 (lower panel). They can hardy 
be detected in the full delay-time signal that is dominated the 
leading lower order moments. I = 0.5x 10 14 W/cm 2 , FWHM 
40 fs, and initial temperature 300 K. 



full period T r . An example containing the effect of the 
fourth cosine-moment is B {td) = ((sin 2 26 1 )) (td), which 
is illustrated in Fig. [7) For a heteromolecular diatomic 
molecule with a permanent dipole moment, the lowest or- 
der interaction Hamiltonian is characterized by the first 
power of cos 9 i.e. n = 1. Thus the alignment measure, 
a cosine moment with N = 2, will show n x N = 2 re- 
vivals within in the full period. Higher order revivals 
may occur since Eq. (|124[) in principle holds for any 
combination (N,n). We may recall, however, that for 
large N, the expectation value might be too weak for the 
lowest fractional revivals to be measured with sufficient 
resolution in practice. This circumstance is illustrated in 
Fig. [5] which shows the high order fractional revivals for 
N = 6 and N — 8 cosine-moments, ((cos 6 #)) (td) and 
((cos 8 #)) (td), along with their magnifications. 



C. Phase Relations of Fractional Revivals 

Can one predict the relative phases of the fractional 
revivals? We may answer this question positively. From 
the phase difference fEq. (|123| )), one finds: 



(Trev) ~ A(j) J 2 %(T rev ) = 4tt 
A0 2 J ° 2 +1 (T rCT /2) - A^ 2 (T rev /2) = 2tt 
A0 2 >O 2 +1 (W4) _ A0 2 }o 2 (T rCT /4) = 7r 



(125) 



Eq. (|125|) predicts that at jT r the phase for J even is an 
exact mirror image of the phase for J dd, as in fact is the 
case in Fig. [51 calculated for N 2 . From the above, we 
may further predict that: 

(a) for O2, which posses J dd levels only, will show a 
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Figure 9: Dependence of alignment moment ((cos 2 9)) (td) 
(for the case of N2) on [Jo, Mo). Panel (a): calculated with 
a single state |5, 0) only (solid line), and retaining all odd Jo 
with J m ax = 19 (dashed line). Panel (b): with |6, 0) only 
(solid line), and all even Jo with J max = 20 (dashed line). 
Panel (c): calculated with all J upto J ma x = 20. Results are 
for pulse intensity 0.8 x 10 14 W/cm 2 and with FWHM 40 fs, 
for initial temperature 300 K. 



levels only, will show a 



"peak" at T r /4, 

(b) CO2, which possesses J e . 
"valley" at T r /4, and 

(c) N2, which possesses both the majority J even levels 
and the and the minority J dd levels in the ratio 2 : 1 
(due to the nuclear statistics of the molecule |22|), will 
show the the revival at T r /4 that would be a "valley" 
like the one for the J eV en levels only, but with only half 
its normal "depth" , due to the counter contribution from 
the minority J D dd levels. 

We note that one may also predict the nuclear statis- 
tics of such molecules by comparing the revival shape 
at T rev /2 and T rev /A. Let us first define a modula- 
tion amplitude at half-revival to be equal to the dif- 
ference between peak and the base (or average) signal: 
(Ax/2 — Sl°2 — S1/2) ■ Similarly, a modulation amplitude 
at quarter revival is equal to the difference between the 
top and the base (average) signal: (A 1 / 4 — S^K — S^Vj). 
The amplitude at half-revival is a sum of even and odd 
J contributions, and therefore A1/2 is always positive. 
In contrast, the amplitude at the quarter-revival arises 
from their difference, and therefore A\u can be positive 
(if it makes a "top" alignment) or negative (if it makes 
an "anti-top" alignment). Therefore, the existence of a 
"top" signal at the quarter-revival is a sign that even J 
levels are dominant. Similarly the presence of an "anti- 
top" signal at the quarter-revival signal is a sign of domi- 
nant odd J levels. From this observation, one can deduce 
the nuclear statistics from the ratio between the effective 
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(finite) number of even and the odd J levels ( J e 
J odd, respectively) excited: 



A 



1/2 



.4 



1/4 



J, 



odd 



A 



1/2 



.4 



1/4 



and 



(126) 



Thus, for example, the dynamic signal of O2 shows 
Al/2 = -^1/4 indicating the absence of the even J lev- 
els. In contrast, Ai/ 2 — —A1/4 for CO2, indicating the 
absence odd J levels. For N 2 , we have Ai/^ = —^Ai/ 2 , 
and hence we have J even ■ J odd = 2:1. This property 
might be used for detecting the existence of isotopes of a 
molecular sample, as has been suggested recently [60]. 



D. Beat Frequencies 



From Eq. (|122|) . it is seen that the phase difference 
associated with (cos 2 9) is (B /h)(AJ +&). For B in cm -1 , 
the phase difference reads 



A</> ( J -> J ± 2) = 2ttBc{AJ + 6) 



(127) 



with c in cm/second. According to Eq. (|127|) , one can 
make a Fourier transform of (cos 2 9) using Be as ba- 
sis frequency and find a series of peaks at (4 J + 6). 
Fig. flTJI shows the Fourier transform of ((cos 2 9)) of N2, 
O2, and CO2. The spectrum of 2 has peak series at 
(10, 18,26, ..)Bc = {AJ dd + 6)Bc, showing that 2 has 
odd J levels only. In contrast, the peak series of CP 
are located at (6, 14, 22, 30, ...)Bc = (4J ei)e „ + 6), sho 
ing that CO2 has even J levels only. For N2, we obta 
a series (6, 14, 22, 30, ...)Bc = (4J eue „ + 6) that is twi 
as strong as the series (10, 18, 26, ..)Bc = (iJ dd + 6)1 
It implies that both even and odd J levels are present 
N2, in the ratio J even : J dd = 2:1. These conclusi 
are consistent with the analysis based on the dynan 
signals. 

For ((sin 2 2$)), there are two kinds of difference 
beat frequency. The first one is related to the transitio 
with A J = ±2 and is expressed by Eq. ([127)) . T 
second one is related to the transitions with A J = : 
and can be expressed as 



A</>( J -> J ± 4) = 2irBc(8J + 20) 



(12 



As a results, in addition to the series of lines (4 J + ( 
the Fourier transform of ((sin 2 20)) also has another s 
ries of lines at (8 J + 20), with A J = 4. Fig. [TTlsho- 
the calculated Fourier transform of ((sin 2 29)) of O2. 
is seen from Fig. [TTJ that the first series (Aj = ± 
reaches its maximum at J max =11, while the second o 
(AJ = ±4) at J m ax = 13. This difference comes frc 
the fact that the AJ = 4 transition requires AJ = 2 
an intermediate transition. As a result, a A J = 4 transi- 
tion can occur one step after the A J = 2 transition; for 
O2 with only J a dd levels present, this implies a a shift in 
J by 2, from J max — 11 to J max = 13, as seen above. 
From Fig. [11] one also finds that the intensity of the sec- 
ond transition is smaller than of the first one. This arises 
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Figure 10: Fourier transform of the alignment moment 
A{td) = ((cos 2 9)) (td), plotted using Be as the basis fre- 
quency (lower scale). Following the (4 J + 6)-rule for the line 
positions from ((cos 2 9)), the peak frequencies are seen to oc- 
cur for odd J only, for the case of O2, even J only for the case 
of CO2, and for both even and odd J, for the case of Na. The 
corresponding J values are shown in upper scale. 
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Figure 11: Fourier transform of ((sin 2 2#)) of O2 : all lines 
(panel a), line from transitions with A J = ±2 only (panel 
b), and lines from transition with A J = ±4 only (panel c); 
pulse of intensity 0.5 x 10 14 W/cm 2 , FWHM = 40 fs, initial 
temperature 300 K. 




3 50 



" 1 1 T | T l TW^T l T | m i T |l |^T | T l TnW l T | T i m i Trffm 

10 20 30 40 50 60 70 80 90 




10 20 30 40 50 60 70 80 90 
Frequency (in units of Be) 

Figure 12: Comparison of the experimental [l3j] vs. the- 
oretical Fourier spectrum for the dynamic 19th harmonic 
signal for N2 (Fig. [4jl. Both spectra show: series 
I: (6,14,22,30,38,..) Be, series II: (10 , 18 , 26 , 34 , 42, ..) Be, 
series III: (20, 28, 36, 44, 52, 60,.. )Bc, and series IV: 
(4,8,12,16,..)Bc. 



from the circumstance that the allowed matrix element 
of the second transition with the greater separation in J 
is weaker than the one with the lesser separation. 



E. Signals in the Frequency Domain 

To further compare with experimental data, we Fourier 
transform the calculated dynamic signals to get their 
spectra in the frequency domain. They may then be com- 
pared with the F.T. of the experimental data. The results 
for the 19th harmonic signal for N2 is compard with the 
experimental data in Fig. [T2j It can be seen that the 
experimental spectrum (panel a) exhibits two prominent 
series I: (6, 14, 22, 30, ..) Be and II: (10, 18, 26, 34 , ..) Be, 
which are also present in the theoretical spectrum (panel 
b). They can be easily understood to arise from the F.T. 
of the ((cos 2 6*)) (td) term in Eq. (|108f> which vanishes 
unless A J = 0, ±2; this produces a sequence of lines 
(Ej +2 - Ej)/2n — (4 J + 6) Be, and gives the series I 
and II, for the even and the odd J levels, respectively. 
The relative prominence of the series I over the series II, 
from both experiment and theory, seen in the two pan- 
els in Fig. Q/21 could be understood as the 2 : 1 ratio of 
the J even over J odd levels, a well-known consequence 
of the nuclear spin statistics of N 2 (e.g. [1, HU). The 
weakly resolved series III: (20, 28, 36, 44, ..)Bc and se- 
ries IV: (4,8, 12, 16, ..)Bc in Fig. \W[a) are the un- 
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Figure 13: Comparison of the experimental []]§] vs. 
theoretical Fourier spectrum of the dynamic 19th har- 
monic signal for O2 (Fig. [5}. Both the spec- 
tra show: series II: (10 , 18 , 26 , 34, 42, ..) Be, series III: 
(20, 28, 36, 44, 52, 60, ..) Be, and series V: (8 , 16 , 24, ..) Be 



expected series that could not be produced by the F.T. 
of the leading term ((cos 2 6*)) (td)- We note that the se- 
ries III and IV, although weak, are certainly also present 
in the theoretical spectrum in Fig. [T2T b). To inter- 
pret their origin, we consider the two higher order terms 

involving ^(cos 2 6>) 2 ^} (td) and ((cos 4 6*)) (td) in the sig- 
nal for N2, Eq. (|108|) . Because of the presence of the 
square of the second moment, the expected beat fre- 
quency from ^(cos 2 #)~^> (td) not only includes the fre- 
quencies (4J + 6)Bc but also their sum and difference 
frequencies, as indicated below: 

(a + b cos u>\t) (a' + b' cos oj[t) 

= aa' + a'bcosuJit + ab' cos ui'yt + bb' cos uj\t cosa^i 

= aa + a'b cos u)\t + ab' cos uj[t 

bb' , bb' . . nn . 

H — — cos (w\ + w 1 ) t + — cos (uii — u-y) t (129) 

Above, the term a arises from transition with A J = 
with u>o — 0. The frequency uj\ arises from transition 
with A J = ±2. The sum frequency (uix + uiy) yields 
(E.j+2 - Ej) /2vr ee (4( J + J') + 12) Be series whereas 
the difference (u>x — lj[) produces (-E/+2 — Ej) /2tt ee 
(4(J - J')Bc) > 0. For integer J and J' they yield 
the series IV: (4, 8, 12, 16, ..)Bc. The next term 
((cos 4 0)) (td) vanishes unless A J = 0, ±2, and ±4 
produces not only (E J+2 — Ej)/2n = (4 J + 6) Be se- 
quences lines but also (-Ej+4 — Ej)/2tt — (8 J + 20)Bc 
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Table III: List of all the predicted series arising from the moments ((cos 2 #)) , <^(cos 2 #) y , and ((cos 4 #)) that are present in 
the expression for the signal for N2 (for which both odd and even J 's are allowed). 



No. Group freq. Weighting factor 


Formula 


Peak series (in Be) 


Expt. series 


«cos 2 *» 


1 


a 









2 wi 


b 


4J + 6 


10,18,26,.... for odd J 
6,14,22,.... for even J 


II 
I 


(<cos^) 2 ) 


3 


ad 









4 uii and uj[ 


a'b and ab' 


4J + 6 


10,18,26,.... for odd J 
6,14,22, for even J 


II 
I 


5 LUl + Uj[ 


bb' 
2 


4(J + J') + 12 


20,28,36,.. 


III 


6 UJ\ — lu[ 


bb' 

2 


4(J — J') > 


4,8, 12,.. 


IV 




7 


a 









8 U)i 


b 


4J + 6 


10,18,26,.... for odd J 
6,14,22, for even J 


II 
I 


9 UJ2 


c 


8J + 20 


28,44,60.. 


III 



Table IV: All possible frequency arising from / (sin 2 cos 2 #) \ for O2 whose only odd J 's are allowed. The weak frequencies 
are noticed with (*). 



No. 


Group freq. 


Weighting factor 


Formula 


Peak series (in Be) 


Expt. series 


1 




aa 









2 


uji and uj'i 


a'b and ab' 


4J + 6 


10,18,26,... 


II 


3 


102 and uj' 2 


a'c and ac' 


8J + 20 


28,44,60,... 


III 


4 


U\ + 0j[ 


bb' 

2 


4(J + J') + 12 


20, 28,36,.. 


III 


5 


UJl — uj[ 


bb' 
2 


4(J — J') > 


8,16, 24,.. 


V 


6 


LU2 + UJ 2 


cc 

~T 


8(J + J') +40 


56, 72,88, .. 


I* 


7 


UJ2 ~ Uj'2 


cc' 


8(J — J') > 


16, 32,48,.. 


V* 


8 


uji + uj'2 and 0J2 + uj'i 


b -f^ b -T 


4(J + 2J') + 26 


38,46,54,.. 


VI* 


9 


UJl — Uj'2 


be' 
2 


4(J-2J')-14>0 


6,14, 22,... 


VI* 


10 


UJ2 — Uj'i 


b'c 

2 


4(-J + 2J') + 14>0 


2,10,18,... 


II* 



gives series III (20, 28, 36, 44, ..)Bc. All the possible 
series arising from these three leading terms and their 
grouping according to those observed experimentally are 
shown in table IIIII Note that series III is identical, 
and overlap, with the series IV: (4 ? 16, ..)Bc and 

adds to its strength. Moreover, the remaining lines at 
(4, 8, 12, 16, 24, 32, ...)Bc found in the experimental 
spectrum in Fig. [T2T a) as well as in the theoretical spec- 
trum in Fig. [T2Tb). confirm the existence of the series IV 
which is distinct from the series III. The existence of se- 
ries III and IV is a prove of the fact that the dynamic sig- 
nal of N2 can not be described in term of ((cos 2 0)) (td) 
only. 

In Fig. [13] we compare the experimental spectrum 
(panel a) for O2 [IH with the theoretical spectrum (panel 
b) calculated from Eq. (| 1 13[) . Both the experimental 
and the theoretical spectra in Fig. [T3] show the Raman- 



allowed series II: (10, 18, 26, 34,42, ..) Be, but not the 
series I: (6, 14, 22, 30, 38, ..) Be, seen for N2. The anoma- 
lous series III: (20, 28, 36, 44,..) Be, discussed in the case 
of N2 above, however, appears for O2 as well. Finally, an- 
other anomalous sequence V: (8 , 16 , 24 , ■•) Be can be 
seen to be present in the data for O2 in Fig. [13k . that, 
we point out, can not be generated by F.T, of ((sin 2 20)) 
term. To interpret the origin of the observed series in O2 
we first consider the leading term given by Eq. (|110p . 

^(sin 2 9 cos 2 #) 2 ^ ■ The matrix element (sin 2 9 cos 2 6) 

vanishes unless A J = 0, ±2, ±4 corresponds to frequency 
ujo, uii, and UJ2 ■ Thus, there will be the various sum and 
difference frequencies that arise from the presence of the 
squared moment, as follows: 



(a + bcosuiit + CC0SCJ2*) ( a ' + b' cosui[t + c' cosoj' 2 t) 
= aa! + ab' cos uj[t + a'b cos ui\t 
+ac cos u>' 2 t + a c cos ioit 



bb' 
T 


cos (wi - 




66' 
2 


cos (wi 






cc' 
2 


COS (W2 - 




cc' 
2 


COS (u>2 


-U)' 2 )t 




6c' 
2 


cos (ui - 




6c' 
2 


COS (wi 






6'c 
2 


COS ()jj 2 ~ 




6'c 

~ T 


COS (W2 


-wi)« 


(130) 



with a > 6 > c. As discussed before, the fre- 
quency wi generates the lines (4J + 6) that for odd 
J give the series II: (10, 18, 26, ..)£?c. The series I: 
(6, 14, 22, ..) Be that would exist for even J is absent 
from the O2 spectrum. This is easily understood 
as due to the nuclear spin of O atoms, which is 0, 
that strictly forbids any even J rotational levels for 
2 (as dictated by the overall symmetry of the to- 
tal wavefunction for O2). For odd J, the frequency 
uj 2 produces the lines (8 J + 20)Bc = (28, 44, 60, ...)Bc, 
whereas the sum frequency (u>i + lj^) produces the lines 
(4(J + J') + 12) Be = (20, 28, 36, ...)Bc\ taken together 
they generate the series III: (20, 28, 36, 44, .,)Bc, Sim- 
ilarly, the difference frequency u>% — uo 2 gives rise to the 
series V: (8 , 16 , 24 , ..)Bc, as shown in table [TVl All 
the above predicted series are observed in the Fourier 
spectrum for O2. It is also shown in table [TVl that fre- 
quency u)\ + lu' 2 and lu 2 + u)[ produces the weak (strength 
of order order "be" ) series VI: (4(J + 2J') + 26)Bc = 
(38, 46, 54, ...)Bc. Despite its weakness the existence of 
this series too is evidenced by the presence of the line at 
38 Be. It is worth noting that the series V and VI can not 
be generated from the moment ((sin 2 2$)) alone. The 
remaining higher order terms in Eq. (|1 13[) contribute, 
generally very weakly, cither to the lines in scries above 
or to some additional lines that can be seen in Fig. [T3T b). 
but hardly resolved in Fig. I13f a). Finally, we may point 
out that the heights of the few lowest frequency lines in 
the experimental data in Fig. [TBTa') for O2 are believed 
to be due to the fluctuation of the laser outputs in the 
experiment (see, foot-note [20] of [II]). 

We may point out that during the test calculations, 
the relative strengths of the lines in a calculated spec- 
trum were found to depend sensitively (cf. Fig. [TJ] ) on 
the assumed molecular temperature, which is rather diffi- 
cult to determine experimentally. This sensitivity, on the 
other hand, provides a way to estimate the temperature 
of the molecular ensemble of interest in the experiment, 
by requiring that the rotational line for the maximum 
height of the Fourier spectrum of the dynamic signal to 
match with the peak of the Boltzmann distribution of 
the initially occupied rotational levels, and adjusting the 
latter to find the matching temperature. 

Fig. [T3] shows a calculated spectrum for initial tem- 
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Figure 14: Calculated spectra for N2 (panel a) and O2 (panel 
b) at a Boltzmann temperature 75 K; laser parameters are as 
in Fig. Hand [12] for N 2 , and, as in Fig. [5]and[rj for Q 2 . 



perature 75 K whose peaks are shifted from one of 200 
K (Figs. [T2l and [T3| . We also point out that our adi- 
abatic theory produces series III in for N2 and series V 
and VI for 02- These terms arise from the cross-term, 
and hence can not be produced from frozen nuclei ap- 
proximation. Thus, the spectrum in frequency domain 
gives more succinct and clearer information of the HHG 
signal and therefore power full to test the model [l9j |. 



Interplay of Polarization Geometry a and Delay 
Time td 



So far we have limited our applications to the HHG 
signal for parallel geometry of the pump and probe po- 
larizations. We now consider the more general case when 
probe polarization is rotated by a given angle a. Fig. 
[TBI (upper panel) shows our computational results of the 
HHG signals as a function of td, at three different fixed 
a, i.e. a = 0°, 45°, and 90°, for N 2 . We note that the 
signal for a = 90° changes its phase by n with respect to 
the signal for a = 0°, a phenomenon that is also observed 
recently fl2l [274. l55l|. In contrast, the signal for a = 45° 
is seen to remain rather flat with change of a. 

To see qualitatively the a dependence of HHG signal 
of N2, we consider the leading term of Eq. (|114p which 
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Figure 15: Calculated 19th harmonic dynamic signal for N2 
(upper panel) and O2 (lower panel) for various pump-probe 
polarization angles, i.e. a — 0°, a — 45°, and a = 90°. The 
laser parameters are similar with one in Figs. [4] and [5] for N2 
and O2, respectively. The initial temperature is 200 K. 



is given by 



— <.W 1 J. n ) 

— c 00 -r c 01 



- sin 2 a + - (3 cos 2 a - l) 
2 2 1 > 



x((cos 2 6)){t d ) 



(131) 



Thus, for 
(t d ;0°) 



the 



parallel polarizations we have, 
°ao + 4i ((cos 2 9)) (i d ) and for 
the perpendicular polarizations, (i<j; 90°) w 

Cqq^ + (l — ((cos 2 0)) which are clearly of op- 
posite phase as a function of t d . These above expressions 
also show that the modulation depth for a = 90° is 
smaller than one for a = 0°, that can not be obtained 
by planar model [12]. Eq. (|13ip also implies that the 
extrema of the signal would occur for sin a cos a = 0, or 
the maximum at a = 0° and the minimum for a — 90°, 
as seen in Fig. I151 upper) and confirmed experimentally 
[12, HE HI- Eq. plTT]) also implies that at a critical 
angle a c given by (3cos 2 a c — l) =0, or a c w 55°, the 
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Figure 16: Variation of the 19th HHG signal as a function of 
pump-probe polarization angle a, near the first half-revival, 
for N2 (upper panel). The observation times are td = 4.090 ps 
for the "top", td = 4.180 ps, the "average", and td = 4.265 ps 
the "anti-top" alignment times. And similarly for O2 (lower 
panel) . The pulse parameters are the same as in Fig. [4] for 
N2 and Fig. [5]for O2. The initial temperature is 200 K. Note 
the existense of a "magic" angle at a = arctan \/2 w 55° for 
N2, where the dynamical signals all coincide, and a "crossing 
neighborhood" near that angle for O2. 



signal essentially remains constant and independent of 
the delay t d between the pulses. This geometry therefore 
can be used to generate a steady state HHG signal from 
N2, with femtosecond pulses. 

The magic angle in fact is a generic signature for the 
a g symmetry of the active molecular orbitals. 

For 2 , the leading term of HHG signal (Eq. (I115j) ) 
reads 



S {n) {t d -a) 



(n) 

-ii- (((-35 cos 4 a + 30 cos 2 a - 3) 
x (cos 4 9) {t d ) 

+ (30 cos 4 a - 24 cos 2 a + 2) (cos 2 9) 

(132) 



the 



parallel 

c 



Thus, for 

4" } ( (sin 2 9 cos 2 9)~J and 
pendicular polarizations, 



3 sin 4 a + 4 sin 2 a) ) \ + 

polarizations we have, 



> ) ((-( C os 4 c9) + (cos 2 e)y 

for the 
S (U; 90°) 



per- 



4^ ((-3 (cos 4 9) +2 (cos 2 6>) + ify It was 

clear that the sign of (cos 4 9) does not change 
and hence the phase of eighth revival also re- 
mains constant, as shown in Fig. 1151 (lower 
panel) and confirmed experimentally [13, E3, f55j | . 
These above expressions also can be expressed as 

S(t d ;90°) w C ^((|(sin 2 c?cos 2 c?)-i(cos 2 fj) + i) 2 ) 

shows that the modulation depth for a = 90° is smaller 
than one for a = 0°, that can not be obtained by planar 
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Figure 17: HHG signals for N2 (upper panel) and O2 (lower 
panel) as simultaneous function of delay time td and the rela- 
tive angle a between the pump and probe polarizations. The 
laser parameters are similar to that in [4] and [5] for N2 and O2, 
respectively. The initial temperature is 200 K. 

model [13 • 

We also calculated the signal at fixed time delay for 
various relative polarization direction between pump and 
probe pulses. In Fig. [16] we show the results for N2 (up- 
per panel) and for O2 (lower panel), near the first re- 
vival td w \T rev . For the case of N2 a coincidence of 
the signals is seen to occur at a critical angle a c w 55°. 
This is fully consistent with the prediction of the near 
id-independence of the signal for N2 at a c , made above. 
The HHG signal (solid curve) at t d = 4.090 ps ('top'- 
alignment) lies above the signal at td = 4.265 ps ('anti- 
top' alignment), for all angles a smaller than the critical 
a c w 55°; the opposite relation holds above a c . For the 
case of O2 (lower panel), in contrast, there is no single 
critical value of a where the signals for all id's could co- 



incide. This is due mainly to the different a-dependence 
of (cos 2 6») and (cos 4 6>) of Eq. (fi~2T3|) . Nevertheless, it 
can be seen that the signal at the 'top' alignment be- 
comes equal to the signal at the 'anti-top' alignment not 
far from a c » 55° (Point a, Fig. [TBT b)). and they reverse 
their relative strengths above it. During the calculation, 
we get that the exact position of points a, b, and c depend 
on the initial temperature. These whole properties well 
agree with the experimental data [26|, [55[ . We also note 
that the calculated a-dependent signal given by Zhou et 
al. gives the same properties for N 2 , but they predicted 
the signal for O2 to be maximized at a w 45°. 

In Fig. [T71 we plot the calculated signals as function of 
both the delay time td and relative polarization angles a 
for both N2 (upper panel) and N2 (lower panel). For both 
molecules, the modulation depths decrease by increasing 
relative polarization angle, reach zero near critical angle 
a c ~ 55° , and increase again but in opposite phase above 
the critical angle. The results for 9 = 90° - 180° are 
exactly mirror image of the results for 9 = 0° — 90° . 

Before concluding this section it is also worthwhile to 
point out that the a-dependence of the HHG signals for 
the more complex tri-atomic molecule CO2 and the or- 
ganic molecule acetylene, HC = CH, because of their ac- 
tive 7T orbital symmetry, are predicted from the general 
structure of the HHG signal given by Eq. (|103p (even 
with out detailed calculations) to exhibit a "cross-over" 
neighborhood near a sa 55°; this is indeed the CcLSC, els has 
been recently observed experimentally [63| . Clearly, the 
presence of the "magic" angle and the cross-over neigh- 
borhood provide a signature of the symmetry of the ac- 
tive molecular orbital, which can be useful in the context 
of the "inverse" problem of molecular imaging [64| from 
the HHG data as suggested first in [20]. Finally, the 
agreement between the present results and experimen- 
tal data provides a clear possibility to control the HHG 
signals by varying both the time- delay and the relative 
pump-probe polarization angle, simultaneously. 



VII. SOME PROBLEMS OF GENERAL 
INTEREST RELATED TO PUMP-PROBE 
SIGNALS FOR HHG 

Before concluding this paper we report on the results of 
our investigations of a number of pump-probe experiment 
related problems of interest in the present context. 



A. Effect of Probe Pulse on the Alignment 

In pump-probe experiments it is generally assumed 
that the dynamical alignment of the molecular axis is 
governed by the ultrashort pump pulse, while the ultra- 
short probe pulse that leads to the HHG signal does dot 
affect the alignment. To check the validity or otherwise of 
this assumption, we directly compare here the dynamic 
alignment moment, A(td',a = 0) calculated as usual as- 
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Delay time (ps) 



Figure 18: Shift of the alignment signal vs. delay time: 
/(cos 2 0)) for N2 (upper panel), and ((sin 2 29)) for O2 (lower 
panel) at 300 K; I pU mp = 0.8 x 10 14 W/cm 2 and I pro be = 
1.7 x 10 14 W/cm 2 with FWHM 40 fs. See, text for further 
explanation. 

suming when only the pump and when both pump and 
probe pulse couple to the molecular polarizability, for N2. 
In the latter case, the total field consists of the superpo- 
sition of the two pulses with a displacement At in time 
between them: 

F(t) = Ft cos (wii) + F 2 cos(w 2 (t- At)) 

+ £20^.92 (t - At) cos (w 2 (t - At)) (133) 

and 

<£ 2 W) = \4o9i(t) + \e 2 w gi{t-At) 
+2£i e 20 5i (t)g 2 (t - At) 
x ((cos (wit)) (cos (w 2 (t - At)))} (134) 

In the above the indices 1 and 2 stand for pump and 
probe pulse, respectively. Suppose the data are recorded 
after the probe pulse dies out, then the observing time is 
t = At + r, where r is the duration of the probe pulse. 
Eq. (fig!)) then reads 

(e 2 (At + T)) = ^e 2 w9l {At + T) + ^e 2 o9l (r) 

+2eio£ 2 o5i (At + t) g 2 (t) 

x ((cos (wi (At + r))) (cos (w 2 r))) (135) 

showing its dependence on the delay between the two 
pulses At and the length of interaction of the probe 
pulse r. In Fig. 1181 we plot the alignment moment 
((cos 2 0» (At + t) and ((sin 2 29)) (At + r) for 2 , plot- 
ted as a function of delay between two pulses At, for a 
fixed r = 40 fs, as shown by the solid curve . The results 



are compared with that obtained from for the pump pulse 
alone (dashed curve), recorded at the same time. The 
comparison clearly shows that the probe pulse changes 
the dynamic alignment ((cos 2 #)) (t^) in that the signal 
is shifted upward by the presence of the probe pulse as 
may be expected from the enhanced intensity of the field 
when both the pulses overlap significantly (before it dies 
out). Thus, except perhaps when the two pulses overlap 
(or are separated only negligibly) this do not change the 
general characteristics of the dynamical signals. There- 
fore, within the above mentioned exception, one may ne- 
glect the effect of the probe pulse on the HHG signal. 



B. Effect of Initial Temperature 

We assume that the rotational eigenstates \JqMq) of 
the molecule are occupied thermally before the interac- 
tion with the pump pulse. Unlike an upward transitions 
(Jo, Mo) — > (Jq,A/o) to the states with an arbitrarily 
high Jq, the downward transition toward Jq > Mq can 
be restricted. As a result, a wavepacket state created by 
the pump pulse would consist of eigenstates with higher 
occupation of Jq > Mo, implying that the vector of ro- 
tational angular momentum would tend to lie in a plane 
perpendicular to the pump polarization direction. Since 
the rotational angular momentum itself is perpendicular 
to the internuclear axis of a linear molecule, the above 
condition, Jq > Mo, means also that the molecular axis 
would tend to align in the direction of the laser polar- 
ization. This is the physical reason why the alignment 
angle of the molecular axis with respect to the polar- 
ization direction after the laser interaction is generally 
smaller after the interaction than before it, i.e. the de- 
gree of alignment increases on interaction with the pump 
pulse. Since at a lower initial temperature, the lower Mo 
states are relatively more occupied initially, the "degree 
of alignment" A = ((cos 2 9)) would tend to be higher, 
allowing the molecules to be more readily aligned at a 
lower initial temperature. 



C. Mean Energy of the Molecule after the Pump 
Pulse 

It is interesting also to examine the way the mean en- 
ergy of the molecule changes with increasing intensity of 
the pump pulse. Fig. [JD] shows the calculated mean en- 
ergy (E)j M (t) at a time t , before and after the arrival 
of the peak of the pump pulse (of length t p ). As expected, 
the figure shows that increasing the peak pulse intensity, 
increases the mean energy of the molecule or, the "effec- 
tive temperature" T e //. = (E} JoMo (t > t p ) /fee, where 
ks is the Boltzmann constant. However, it should be re- 
membered that after the pulse interaction, the molecular 
system is not a state of thermal equilibrium, rather it is 
in a state of dynamical equilibrium (or steady state) that 
can not be characterized thermodynamically. 
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Figure 19: Influence of the initial temperature, 50° or 300°, 
on the molecular axis-distribution (panel a-b), on the lead- 
ing dynamic moment ((sin 2 2(9)) (td) (panel c-d), and the 
Fourier spectrum (panel e-f) of the latter, for the case of 
O2. The molecular axis-distribution (panel a-b) are given 
at three values of the delay time td, i.e. at the "top" align- 
ment td ~ 5.649 ps (dashed line), at the "average" alignment 
td = 5.812 ps (solid line), and at the "anti-top" alignment 
td = 5.975 ps (dotted line); the radii are in the same scale. 1= 
0.5 x 10 14 W/cm 2 and FWHM=40fs. 



To estimate an effective "temperature" of the rota- 
tional wavepacket states \§j M (t)}, in the steady state 
regime, i.e. for (t pump <t< t pro h e ), we note (a) that 
the rotational wavepacket states \<&j Q M {t)) form a lin- 
early independent set of states like the set of rotational 
eigenstates | JqMq) from which they evolve, (b) that the 
individual rotational wavepacket states evolve in one-to- 
one correspondence with the initially occupied rotational 
eigenstates | JqMq), (c) that the mean energy of each of 
the rotational wavepacket states reach a steady state, 
also one-to-one of energy from (E) JqMq (t > t pump ) = 
Ej m + (E) JoMo (t> t pum p). If further the above 
change in the mean energy (E) JoMq (t > t pump ) - 
(E) j M (t > tp Ump ) is independent of the the individual 
states chosen (indices { Jo-Mo}) then one might use it to 
define an effective "temperature" change, AT e ff, given 
by 

= (AE) (t>t pump ) 

KB 

We may note in Fig. [23 that the change in the mean 
energy in the steady state regime is indeed essentially 
independent of the states of the system chosen. Thus, 
the effective "temperature" of the system, at the end of 
the interaction with the pump pulse, becomes 

T eff = T + AT e ff (137) 

Note that T e f /. is in general greater than the initial gas 
(jet) temperature, To, and it tends to increase with the 
increase of the pump intensity. Later on we shall describe 
a method of determining this "effective temperature" of 
the system from a theoretical analysis of the experimental 
HHG data. 



D. Some Non-equivalent Definitions of the HHG 
Signal 
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Figure 20: Mean energy of O2 before and after interaction 
with the laser pulse, for different initial temperatures. The 
vertical dashed lines indicate the extent of the pulse duration; 



1= 0.5 x 10 J 
explanation. 



W/cm 2 FWHM=40fs. See, text for further 



In this sub-section we briefly discuss two alternative 
definitions of HHG signals that have been employed 
earlier and compare them with the definition of the 
HHG signal of the present theory, and with experimen- 
tal data. The present theory defines the quantum tran- 
sition amplitudes for the linearly independent reference 
states \xi(t) >,i = {e, JqM }, (consisting of the prod- 
uct of the ground electronic and the coherent rotational 
wavepacket states) to obtain the independent harmonic 
emission probabilities, and in accordance with the quan- 
tum statistical theory averages the latter to define the 
HHG signal (cf. e.g. Eq.(62]): 



(t d ) = 

CJ2j oMo P(Jo) (*JoMo{U,0) Ti n) {6) <s>j oMo (t d ,e) 



(138) 

It is worth noting that the quantum amplitude calcula- 
tion in the present theory corresponds to the "adiabatic 
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Figure 21: Theoretical Fourier spectrum of the dynamic 19th harmonic signal for O2; pump intensity I = 0.5 x 10 14 W/cm 2 , 
probe intensity / = 1.2 x 10 14 W/cm 2 W/cm 2 , duration 40 fs, wavelength 800 nm, and temperature 200 K. The calculations are 
done using the present theory, Eq. (|138p . (left-lower panel), model A Eq. (|139p . (right-upper panel), and model B, Eq. (|140[) . 
(right-lower panel). For comparison, the experimental spectrum (left-upper panel) is also shown. 



nuclei" approximation [611 162| . in which the matrix el- 
ements with respect to the rotational wavepacket states 
are evaluated at the level of the adiabatic amplitude- 
operator, T^ n \ff), and not at the level of the adiabatic 

probability- operator, |T(")(e)| 2 , that occurs in the more 
drastic "frozen nuclei" approximation. In this theory, as 
in the laboratory, the operational angle is the relative 
polarization angle a, and not the angle between the po- 
larization direction and the molecular axis, 0. In fact, 
the angle 9 is a coordinate that is, as appropriate for a 
quantum formulation, to be integrated over to obtain the 
quantum transition amplitude with respect to the rota- 
tional wavepacket states. 

In the present notation, the two other definitions of the 
HHG signal that have been used earlier (to be referred 
to below as A and B) are, (i) definition A (cf. Eq. (22) 
of [11 and Eq. (6) of ||): 



my 



Jo M 



9) 



4 n) fa) = 

J M (td, I 

Jo M 



rj n )(0,o) 



and, (ii) definition B (cf. Eq. (12) of \Z§ and Eq. (4) of 



Clearly the HHG signals according to models A and B 
differ with each other, and they differ from the present 
definition, Eq. (|138p above. 

We note that model A, Eq. (|139[) . defines the statis- 
tically averaged signal by weighting the individual am- 
plitudes (!) first, and then taking the absolute square 
of the weighted sum. This runs counter to the quantum 
statistical theoretical approach of averaging the probabil- 
ities (not amplitudes) and/or the expectation values of 
Hermitian obscrvablcs themselves, and not their Fourier 
transforms (that are proportional to the emission am- 
plitudes). Furthermore, the definition of model A (Eq. 
(|f 39p ) makes the signal to depend on the mixed products 
of the statistical weights that are in principle indepen- 
dent. 

The signal defined by model B, Eq. (I140j) . is seen to 
depend on the weighted sum of the diagonal matrix ele- 
ments (between the rotational wavepacket states) of the 

" probability operator" |t(™)(#)| 2 - this, of course, is 
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Figure 22: Calculated 9th HHG spectrum of N2 for various 
pump-probe polarizations angle a — 0° and a — 90°; pump 
intensity I = 2x 10 12 W/cm 2 , duration 300 ps; probe intensity, 
I = 5 x 10 14 W/cm 2 , duration 70 fs, and wavelength 798 nm; 
Boltzmann temperature 25 K. 



pump pulse and a short (70 fs) probe pulse, as in an 
adiabatic alignment experiment [6 71 ] used earlier. The 
results of our calculations for N 2 , using Eqs. (|114[) for 
both a = 0° (solid curve), and a = 90° (dash-dot curve), 
are shown in Fig. [22] For the sake of comparison we also 
show the intensity profile (dashed curve) of the pump 
pulse (right scale). As can be seen immediately from the 
figure, in the parallel case, the HHG signal closely follows 
the evolution of the long pump pulse itself (which might 
be expected for an adiabatic process) and the maximum 
of the signal occurs at the maximum of pulse, for a = 0°. 
On the other hand, a minimum is predicted for the signal 
at the maximum intensity, in the perpendicular case, for 
a = 90°. These characteristics of the adiabatic signals 
for N2 are consistent with the experimental observations 
made some time ago [13, IH, . 



not equal to the weighted sum of the absolute squares of 
the diagonal matrix elements of the transition operator 
T^ n \9). In model B the above circumstance is a conse- 
quence of the more drastic "frozen nuclei" approximation 
and an effective inclusion of all transitions, those between 
the same wavepacket states ("elastic-like") as well as 
those between the different wavepacket states ( "inelastic- 
like"). However, unlike in model A (Eq. (|139p ). in model 
B (Eq. (I140|) l the weighted statistical sum is taken, in 
accordance with the quantum statistical theory, at the 
level of the probabilities. 

In Fig. [5T]we show a comparison of the experimental 
data (upper-left panel) of the F.T. of the dynamic signal 
(q = 0) of the 19th harmonic, for the case of O2 with 
the results of the test calculations from: the present the- 
ory theory (lower- left panel) using Eq. (|138[) . model A 
(upper-right panel) using Eq. (|139p . and model B (lower- 
right panel) using Eq. (|140[) . keeping everything else the 
same - the parameters are the same as in Fig. 1131 The 
similarity of the calculated spectrum from the present 
theory and the experimental data is seen to be very sat- 
isfactory, that from model A is similar, except that the 
ratio of the peaks of the series III to series II is gener- 
ally too small compared to the experimental data, and 
series V (8 , 16 , 24 , ••) Be is rather weakly developed. 
In the case of model B, on the other hand, the series 
V (8, 16, 24, ••) Be is simply missing. It may be noted 
that the present comparison also illustrates the ability 
of the experimental data at the level of the (discrete) 
Fourier spectrum, to better distinguish between the var- 
ious theoretical models than may be possible at the level 
of the time-domain signal. 



E. Adiabatic Alignment 

Finally, we apply the present dynamic theory also to 
the adiabatic case, in which we choose a long (300 ps) 



VIII. CONCLUSIONS 

To conclude, we have presented an ab initio intense- 
field S-matrix theory of dynamic alignment of linear 
molecules and the characteristic HHG signals from them 
as detected in intense-field femtosecond pump-probe 
experiments. Useful analytical expressions for the molec- 
ular alignment and the HHG signal as a function of both 
the delay-time, td, and the relative polarization angle, a, 
between the pump and probe pulse, are derived. Thus, 
we give the general HHG signal Eq. (|103p . the signal for 
N2 (generically, active a g orbital symmetry) Eq. (|114[) . 
or its leading term Eq (|13ip . as well as the signal for 
O2 (generically, active ir g orbital symmetry) Eq. (|115p . 
or its leading term, Eq. (|132p . They are used to make 
detailed analysis of the molecular alignment and the 
observed experimental data for the HHG signals from 
coherently rotating N 2 and O2 molecules, both in the 
time-domain and in the frequency domain. The results 
show a remarkable agreement between the theory and 
the experimental observations. Additional predictions 
about the existence of critical relative polarization 
angles, a Cl and their relation to the symmetry of the 
active orbitals and the form of the dynamic signals are 
made. At a "magic" angle, a c « 55°, the dynamic HHG 
signals for all delay times td, are predicted to approach 
each other closely for a linear molecule with a a g orbital 
symmetry, or exhibit a "crossing neighborhood" , for 
the -Kg orbital symmetry; it is also predicted to produce 
a steady emission of high harmonic radiation at the 
magic angle from N2. Moreover, we have investigated a 
number of theoretical questions and experimental effects 
of general interest in connection with the interpretation 
of the pump-probe HHG signals. Finally, we have shown 
that the case of "adiabatic-alignment" , and the resulting 
HHG signal, can be analyzed and understood equally 
well within the present dynamical theory, using simply 
a long duration of the pump pulse. 
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